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I.   INTRODUCTION 


A  model  of  gas  exchange  in  the  avian  lung  is  developed.   The  physical 
structure  of  the  model  is  based  upon  that  proposed  by  Zeuthen  [1].   Dif- 
ferential equations  are  derived  to  describe  the  function  of  the  model.   In 
addition  to  gas  transfer  between  blood  and  air  which  occurs  across  the 
exchange  surfaces  of  the  lung  and  gas  movement  by  bulk  flow  through  the  lung, 
the  mathematical  description  of  the  function  of  the  model  also  accounts  for 
movement  of  a  gas  component  through  the  air  spaces  of  the  lung  by  diffusion. 
Cases  which  treat  diffusion  in  the  air  capillaries,  diffusion  in  the  para- 
bronchial  lumen,  and  diffusion  in  both  these  spaces  are  considered.   In  those 
cases  which  describe  diffusion  through  the  air  spaces,  the  dissociation  curve 
relating  concentration  and  partial  pressure  in  the  blood  is  linear.   Another 
case  has  been  considered  in  which  the  diffusion  through  the  air  spaces  is 
neglected,  but  a  more  accurate  representation  of  the  non-linear  oxygen 
dissociation  curve  is  made.   By  approximating  that  curve  as  a  polynomial  of 
second  or  third  order,  integration  of  the  differential  equations  describing 
gas  transfer  between  air  and  blood  capillaries  can  still  be  accomplished. 
This  latter  technique,  although  applied  here  to  the  avian  lung,  might  also 
find  application  in  analysis  of  gas  transfer  between  pulmonary  capillaries  and 
the  alveoli  of  mammalian  lungs.   With  diffusion  not  accounted  for,  and  with 
a  linear  dissociation  curve,  the  model  can  be  shown  to  be  equivalent  to 
that  of  Scheid  and  Piiper  [2]. 

Solutions  of  the  differential  equations  are  evaluated  and  the  results 
displayed  graphically.   Under  different  conditions  of  cardiac  output  and 
ventilation  the  partial  pressure  profiles  of  oxygen  and  carbon  dioxide  in 
the  blood  capillaries,  air  capillaries,  and  parabronchial  lumen  are 


displayed.  The  cases  of  diffusion  and  no  diffusion  are  contrasted.  The 
case  of  the  model  with  a  non-linear  dissociation  curve  has  also  been  utilized 
to  briefly  analyze  performance  of  oxygen  uptake  by  the  lung. 

The  effort  which  resulted  in  the  avian  lung  model  presented  in  this 
thesis  originally  arose  out  of  a  speculation  on  how  to  determine  possible 
locations  for  the  C0~  receptors  (see  reference  [3]  for  a  description  of 
these  receptors)  which  are  present  in  the  avian  lung.   Single  unit  recordings 
from  the  vagus  nerve  yield  information  about  the  CO.  partial  pressure  in  the 
microenvironment  of  the  receptors  themselves.   If  a  correspondence  could 
be  established  between  CO  partial  pressures  and  various  locations  within 
the  lung  under  given  conditions  of  cardiac  output,  ventilation,  and  partial 
pressure  in  the  venous  blood  and  inspired  gas,  then  the  possible  locations 
could  be  mapped  out.   Considerable  economy  of  effort  would  be  achieved  by 
confining  the  search  for  these  receptors  to  the  indicated  regions. 

The  function  of  the  model  is  based  upon  a  parabronchial  structure  which 
contains  a  representation  of  air  capillaries  and  blood  capillaries  in 
addition  to  the  lumen,  since  the  receptors  might  be  located  within  any  of 
these  structures.   Since  the  physical  dimensions  of  air  capillaries,  blood 
capillaries,  and  parabronchial  lumina  are  so  small  that  experimental 
measurements  of  partial  pressure  cannot  be  made  within  them,  a  mathematical 
model  provides  a  method  of  determining  their  partial  pressures. 

The  actual  determination  of  possible  locations  for  the  C0_  receptors 
has  not  been  performed,  since  the  experimental  results  which  have  been 
published  are  not  in  a  form  which  allows  comparison  with-  the  results  of 
the  model. 

The  form  in  which  the  model  was  initially  developed  is  that  presented 
in  Cases  1,  2,  3,  4  of  Chapter  IV.   As  computations  were  made  it  became 


apparent  that  the  assumption  of  a  linear  dissociation  curve,  although  valid 
for  carbon  dioxide,  was  resulting  in  unrealistic  results  for  oxygen  exchange. 
So  a  secondary  objective  became  the  development  of  a  simple,  yet  reasonably 
accurate,  way  of  incorporating  the  sigmoid  dissociation  curve  of  oxygen  into 
the  model.   This  objective  was  met.   The  solution  is  Case  5  of  Chapter  IV. 


II.   THE  MODEL 

The  mathematical  model  which  is  developed  here  is  fundamentally  a  model 
of  the  avian  parabronchus.   It  is  assumed  that  several  hundred  of  these 
parabronchi,  all  identical,  function  "in  parallel"  to  make  up  the  avian  lung. 
These  parabronchi  are  all  to  be  of  identical  structure,  to  receive  identical 
flows  of  blood  and  gas  which  enter  the  parabronchi  at  identical  compositions. 
This  is  an  accurate  representation  of  the  paleopulmo,  but  not  of  the  neopulmo, 
which  is  composed  of  parabronchi  of  widely  varying  lengths.   But  since  the 
neopulmo,  when  present,  makes  up  a  small  fraction  of  the  lung,  exceeding  25% 
in  no  species,  it  is  perhaps  not  unreasonable  to  assume  the  paleopulmonic 
construction  for  the  whole  lung.   Assuming  no  blood  is  shunted  past  the  lung, 
the  rate  of  blood  flow  to  a  single  parabronchus  is  given  by  the  cardiac 
output  divided  by  the  number  of  parabronchi.   The  ventilation  of  an  individual 
parabronchus  is  assumed  to  be  given  by  the  minute  ventilation  divided  by 
the  number  of  parabronchi.   This  implies  that  all  the  gas  included  in  the 
figure  for  minute  ventilation  makes  one  pass  through  a  parabronchus. 

In  a  bird  the  parabronchus  is  basically  a  cylindrical  structure  which 
is  penetrated  by  a  lumen  along  its  central  axis.   Ventilated  gas  passes  down 
the  lumen.   Along  the  wall  of  the  lumen  are  many  small,  transverse  bands  of 
smooth  muscle  which  define  the  walls  of  air  pockets  termed  "atria".   The 
atria  open  into  the  lumen,  and  evaginating  radially  outward  from  them  are 
many  air  capillaries.   Coursing  adjacent  to  the  air  capillaries  are  the 
blood  capillaries.   The  direction  of  blood  flow  in  these  is  from  the 
arterioles  at  the  outer  periphery  of  the  parabronchus  inward  toward  the  venules 
located  adjacent  to  the  lumen.   Surrounding  the  parabronchial  lumen,  this 
network  of  air  and  blood  capillaries  together  makes  up  the  parabronchial 


mantle.   It  is  in  the  mantle,  across  the  surfaces  separating  blood  from  gas 
that  oxygen  and  carbon  dioxide  exchange  takes  place.   For  a  detailed 
description  of  avian  lung  anatomy,  see  Dunker  [4,5]. 

The  physical  structure  of  the  model  which  is  to  represent  the  actual 
parabronchus  is  illustrated  in  Figure  1A.   This  representation  in  essence 
is  due  originally  to  Zeuten  [1].   In  it  the  detailed  structure  of  smooth  muscle 
bands  and  atria  are  ignored,  and  the  air  capillaries  are  assumed  to  radiate 
outward  directly  from  the  parabronchial  lumen. 

The  function  of  the  model  reflects  the  function  of  a  real  parabronchus. 
Venous  blood,  with  a  high  concentration  of  CO.  and  a  low  concentration  of 
02,  departs  from  a  large  vessel  through  a  series  of  blood  capillaries.   The 
venous  blood  enters  all  blood  capillaries  with  identical  composition.  As 
each  blood  capillary  runs  alongside  its  corresponding  air  capillary,  CO, 
diffuses  from  the  blood  capillary  into  the  air  capillary  and  0.  diffuses 
from  the  air  capillary  into  the  blood  capillary.  After  undergoing  gas 
exchange,  the  blood  from  all  blood  capillaries  flows  into  a  large  vessel  and 
mixes  there  to  form  arterial  blood.  Within  the  air  capillary  there  is  no 
flow  of  gas;  mass  transport  occurs  by  diffusion.   The  diffusion  of  gases 
occurs  so  that  the  0^   supply  is  constantly  replenished  from  the  parabronchial 
air  and  CO^  constantly  drains  into  the  parabronchial  air.  Therefore,  as  gas 
flows  down  the  parabronchial  tube  it  loses  0„  and  picks  up  CO,. 

Explicitly,  the  more  important  assumptions  of  the  model  are: 

1.  All  parabronchi  composing  the  lung  are  structurally  and  functionally 
identical. 

2.  The  model  operates  in  steady  state.   Transient  changes  in  ventilation, 
such  as  those  occurring  during  the  respiratory  cycle,  cannot  be 
handled.  Nor  can  transient  changes  in  the  partial  pressure  of 

gas  or  blood  entering  the  lung  be  handled. 


VENOUS  BLOOD 


compartment  n 


Figure  1.   A:  The  parabronchial  model,  modified  after  Zeuthen.   B:  Detail  of  air 
capillary  and  blood  capillary  compartments.   Expressions  describing  transport  of 
the  considered  gas  between  compartments  are  shown.   C:  Detail  of  parabronchial 
luminal  compartments.   Expressions  are  shown  which  describe  transport  between 
luminal  compartments  by  both  bulk  flow  and  diffusion  and  transport  between  a 
luminal  compartment  and  the  air  capillaries  which  open  into  it. 


3.  The  entire  length-  of  the  parabronchus  is  equally  perfused  with 
blood. 

4.  The  gas  exchange  surface  is  homogeneous.   That  is,  it  has  equal 
diffusion  capacity  per  unit  area  everywhere. 

5.  The  physical  structure  of  the  atria  and  bands  of  smooth  muscle 
is  ignored. 

6.  No  blood  is  shunted  past  the  lungs  without  exposure  to  gas  exchange 
surfaces. 

7.  There  is  no  interaction  between  gas  components.   The  Bohr  and 
Haldane  effects  are  not  accounted  for. 

The  model  developed  in  this  thesis  generalizes  in  several  important 
ways  upon  the  model  developed  by  Scheid  and  Piiper  [2].   In  their  model, 
Scheid  and  Piiper  neglected  the  structure  of  the  air  and  blood  capillaries. 
Gas  exchange  is  represented  as  taking  place  between  a  thin  sheet  of  blood 
flowing  across  the  surface  of  a  tubular  parabronchus  and  the  gas  flowing 
through  that  tube.   The  direction  of  blood  flow  is  perpendicular  to  the 
direction  of  gas  flow,  giving  rise  to  the  description  of  their  model  as 
"cross-current".   In  the  model  developed  in  the  following  pages,  the  air 
and  blood  capillaries  are  retained  as  separate  structures,  and  expressions 
are  developed  which  describe  the  variation  in  partial  pressure  of  a 
particular  gas  component  along  the  length  of  both  the  blood  capillaries  and 
the  air  capillaries.   Also,  in  the  model  of  Scheid  and  Piiper  gas  is  assumed 
to  move  through  the  parabronchus  by  bulk  flow  alone,  without  diffusion  of 
any  gas  component  in  the  axial  direction;  in  the  model  developed  here  both 
bulk  flow  and  diffusion  are  accounted  for.   A  further  generalization  in  the 
developed  model  lies  in  the  handling  of  the  non-linear  oxygen  dissociation 


curve.   In  Che  model  of  Scheid  and  Piiper,  the  oxygen  dissociation  curve 
has  been  represented  as  linear  in  order  that  the  differential  equations  may 
easily  be  integrated.   However,  it  has  been  found  here  that  by  handling 
those  equations  in  a  slightly  different  way  and  representing  the  dissociation 
curve  by  a  polynomial  of  second  or  third  order,  the  integration  may  still 
be  performed.   Use  of  a  polynomial  allows  a  much  more  accurate  representation 
of  the  dissociation  curve  over  the  physiological  range. 


III.   DERIVATION  OF  DIFFERENTIAL  EQUATIONS 

3.1  The  Air  Capillary-Blood  Capillary  Pair 

The  diffusion  occurring  between  an  air  capillary  and  a  blood  capillary 
is  considered.   Although  a  mixture  of  different  gases  is  present,  it  will 
be  assumed  that  they  do  not  interact.  Attention  can  therefore  be  focused 
on  a  single  component  of  the  mixture. 

Both  the  air  capillary  and  blood  capillary  are  represented  as  a  string 
of  homogeneous  cells  with  mass  transfer  occurring  across  cell  surfaces. 
Figure  IB  illustrates  this  geometry. 

A  blood  capillary  cell,  numbered  n,  will  undergo  mass  transfer  with 

blood  capillary  cells  n-1  and  n+1,  and  with  air  capillary  cell  n.   The  rate 

at  which  the  gas  component  under  consideration  diffuses  into  blood  capillary 

cell  n  from  air  capillary  cell  n  is 

hAx(P    -  P.   )   •  (3.1-1) 

g,n    b,n 

The  rate  at  which  this  component  is  transferred  due  to  flow  from  blood 

capillary  cell  n  to  blood  capillary  cell  n+1  minus  the  rate  at  which  the 

component  is  transferred  from  blood  capillary  cell  n-1  to  blood  capillary 

cell  n  is 

f(Cb,n  "  Cb,n-1>  ■  t3-1"2) 

Under  steady  state  conditions  (3.1-1)  and  (3.1-2)  are  equal: 

f(Cb,n  "  S.n-l5  "  h^g,n  "  Pb,n>Ax  '         ".1-3) 
Dividing  both  sides  of  (3.1-3)  by  Ax,  then  taking  the  limit  as  Ax-K)  leads  to 

3T  +  K"K  •  ".1-4) 
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In  Cases  1,  2,  3,  and  4  the  simplifying  assumption  is  made  that  C,  and 

b 


(3.1-5) 

:  P,  or  i 

D 

Equation  (3.1-5)  substituted  into  (3.1-4)  yields  the  differential  equation 

dlT  +  If  Pb  =  W  Pg   •  "-1-6) 

Since  it  is  assumed  that  the  state  of  venous  blood  entering  the  lung  is 
known,  the  initial  condition  for  (3.1-6)  is 

Pb^0)  =  Pb,v   '  °r    Cb^=0)  =  S,v   •       (3-W) 

An  air  capillary  cell,  numbered  n,  will  undergo  mass  transfer  with  air 

capillary  cells  n-1  and  n+1,  and  with  blood  capillary  cell  n.  The  rate  at 

which  the  component  diffuses  into  air  capillary  cell  n  from  blood  capillary 

cell  n  is 

hAx(P,    -  P   )  •  (3.1-8) 

b,n    g,n 

In  the  actual  air  capillary  where  concentration  of  any  gas  component  is  a 

continuously  changing  quantity,  the  rate  at  which  a  component  diffuses  across 

a  point  in  the  air  capillary  is 

dP 
-cn^*  ■  (3.1-9) 

In  Che  discrete  cell  model,  the  rate  at  which  the  component  diffuses  from 
cell  n  to  cell  n+1  is 

-gri  ?s,n+L" ?s,n  ■  "-1-10) 

The  rate  at  which  the  component  diffuses  from  cell  n  to  cell  n-1  is 
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-«n  ?8'°-^"  ?^n  .  (3.1-11) 

Under  steady  state  conditions,  (3.1-8)  can  be  equated  to  the  sum  of  (3.1-10) 
and  (3.1-11).   Upon  rearranging  this  becomes 


ax  [ 


g,n+l  "  g,n  _  s,n  ~  g,n-lj  _  h_      =  _  h_ 

Ax  Ax      J   an.  g,n     on  b,n 


(3.1-12) 
Taking  the  limit  Ax  +  0, 

*\       h 

For  the  problem  of  diffusion  between  an  air  capillary  and  a  blood  capillary, 
it  is  assumed  that  the  state  of  the  gas  at  the  air  capillary  opening  is 
known.   It  is  also  assumed  that  the  concentration  gradient  at  the  closed 
end  of  the  capillary  is  zero,  since  there  is  no  diffusion  across  the  closed 
end.  Hence  the  boundary  conditions  for  (3.1-13)  are: 

Pg(X=£)  =  Pg,p  =  Pp  '  (3.1-14) 

dP  (x=0) 
aad  —^ 0  .  (3.1-15) 

A  single  compartment  composed  of  a  segment  of  an  air  capillary — blood 
capillary  pair,  as  shown  in  Figure  2,  is  considered. 

The  rate  at  which  the  gas  component  enters  the  compartment  at  the 
upstream  point  is 

dP 
Bf  Pb,u-arl-dx^  •  O-1"16) 

The  rate  at  which  the  gas  component  leaves  the  compartment  at  the  downstream 

point  is 


-»► 

BLOOD    CAPILLARY                               |       "^ 

^      AIR  CAPILLARY                                  |       ** 

1 

Figure  2.   Single  compartment  composed  of  a  segment  of  a  blood- 
air  capillary  pair.   The  upstream  and  downstream  ends  of  the 
compartment  are  at  x„  and  x^,  respectively. 
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dP 

8f  pb,d  - a*  ^r  •  "-1-17) 

Under  steady  state  conditions  (3.1-16)  equals  (3.1-17): 

fdP  dP       i 

Bf (pb,d  -  pb,u>  ■  *«[-$r  -  -fcrj  •         (3-1-18) 

Now  letting  the  upstream  point  be  at  x=0,  and  the  downstream  point  be  at  x=l 
causes  (3.1-18)  to  become 

dP 

Bf  (P.  .  -  P.  J  =  an  —&*£■     ,  C3.1-19) 

b,a    b,v        dx 

where  the  boundary  condition  (3.1-15)  has  been  used.   Thus,  the  rate  at  which 

a  single  air  capillary  leaks  the  gas  component  into  the  parabronchus  can  be 

expressed  either  as 

Sf  (Pb,v  "  Pb,a>  '  °-1-20) 


dP 

-J*  •  (3.1-21) 


3.2  The  Parabronchus 

As  gas  flows  through  the  parabronchus  the  concentration  of  the  gas 
component  changes  due  to  mass  transfer  between  the  parabronchus  and  the 
air  capillaries.  The  model  geometry  to  illustrate  this  process  is  shown 
in  Figure  1C. 

By  a  compartmentalization  similar  to  that  used  on  the  capillaries, 
the  parabronchus  is  considered  to  he  divided  into  a  string  of  cells  which 
are  numbered:   — ,n-l,n,n+l, —  They  are  all  to  have  thickness  iz. 

The  rate  at  which  the  gas  component  is  transferred  due  to  flow  from 
parabronchial  cell  n  to  cell  n+1  minus  the  rate  at  which  the  component  is 
transferred  due  to  flow  from  cell  n-1  to  cell  n  is 
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F(C    -  C    ,)  .  (3.2-1) 

p,n    p,n-l 


This  can  be  rewritten  as 


BF(P    -  P    ,)  .  (3.2-2) 

p,n    p,n-l 

The  rate  at  which  the  component  diffuses  from  parabronchial  cell  n  to  cell 
n+1  plus  the  rate  at  which  it  diffuses  from  cell  n  to  cell  n-1  is 

"*8  Pp,afte"  ?P,n  +  *8  ^^  a/"'1  "         (3.2-3) 
The  rate  at  which  the  gas  component  is  transferred  to  the  parabronchial  cell 
from  the  air  capillaries  is  found  either  from  (3.1-20)  to  be 

wSf(Pb,v  "  Pb,a)Az    '  (3-2_4) 

or  from  (3.1-21)  to  be 

dP 
wan  — jji£  Az   •  (3.2-5) 

Under  steady  state  conditions  (3.2-4)  can  be  equated  to  the  sum  of  (3.2-2) 

and  (3.2-3).   Then  dividing  by  Az  yields 

P..  -P      P    -P    ,-,      ,-P    -P 
P."-+1 EaE  _  P."    P. n-1   ,   r  p,n    p,n-l 


Az  {  )+   BFr~   Az'    j  "  w8f(Pb,v  "  Pb,a) 

(3.2-6) 
Taking   the  limit  Az-K)  yields   the  differential   equation 

9  ^ -¥■£■•  *¥■<?*      "ph     >    •  (3.2-7) 

,2  <p      dz  9  b,a  b,v 

Similarly,  equating  (3.2-5)  to  the  sum  of  (3.2-2)  and  (3.2-3)  leads  to  the 
differential  equation 

A2?  dP  dP 

9_Ji  _  BF  _p_   _  war,  _£ip_ 

,2    *  dz       t>  dx  (3.2-8) 
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The  boundary  conditions  on  (3.2-7)  and  (3.2-8)  are 

pp(z=0)  =  Pp,i  •  (3.2-9) 

dP  (z=X) 
and  ~SJi =0  '  (3.2-10) 

A  single  compartment  defined  by  the  region  between  two  cross-sectional 
cuts  through  the  parabronchus  perpendicular  to  its  axis  is  considered.  The 
compartment  will  contain  the  parabronchial  lumen  together  with  segments 
of  the  venous  and  arterial  blood  vessels  as  shown  in  Figure  3.   The  rate 
at  which  the  gas  component  enters  the  parabronchial  lumen  at  the  upstream 
location  is 

dP 

FBPP,u  - 9*  -if  ■  <3-2-11) 

The  rate  at  which  it  leaves  the  parabronchial  lumen  at  the  downstream 

location  is 

dP] 
'p,d   "T  ~~ dz 

The  rate  at  which  the  dissolved  gas  component  enters  the  compartment  by 

way  of  arterial  and  venous  vessels  is 

S(i,A,v  +  Wb.a.d)   •  «-2-13> 

The  rate  at  which  the  dissolved  gas  component  leaves  the  compartment  by 
way  of  blood  vessels  is 

6(5   P,    +  Q   P,     )   •  (3.2-14) 

^v,u  b,v   Ha,u  h,a,u  *•      ' 

Under  steady  state  conditions  the  quantity  of  gas  entering  the  compartment 

equals  the  quantity  of  gas  leaving  the  compartment,  which  leads  to 


'BPn .H  "  9*  ~%r     ■  (3.2-12) 
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dP 

FBP    -  9<i>  — 1^  +  8(0   .P.    +  Q   .P.    ,) 
p,u        dz      T?,d  b,v    a,d  b,a,d 

dP 


-  FBP    -  8*  — |^  +  8(0,   P.    +  Q   P.     )  • 
p,d        dz      T?,u  b,v    a,v  b,a,u 

(3.2-15) 


P    -P.i 

p,u    p.i 

P    =  P 
P>°    p,e 

dP      dP   . 
P>u  =   P.1   , 
dz      dz 

dp  A 
dz 

Vu=Vd=Q   ' 

P      ■  P 
b , 2 , u    b,a 

Substitution  of  these  boundary  conditions  into  (3.2-15)  and  rearrangement 

yields  the  following  expression  for  P 

For  convenient  reference  the  differential  equations  and  boundary 
conditions  are  listed  in  Table  3-1. 


Table  3-1 
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-^  +  ^P  =  *P 
dx    f  b   f  g 


^  +  ^P  .h,, 

dx    gf  b   gf  g 


Cb(x=0)    - 

:S,v 

P,  (x=0)   = 

0 

:Pb,v 

dx2     "°r 

-  P    ■  - 

I      g 

CtH-       D 

P    (x=Jt)  .- 

i  P 
g.P 

p 
p 

dP    (x=0) 

8 

=  0 

dx 

d2P 

$     dz 

<j>      v  b,a 

"  Vv> 

d2P 

m,    dP. 

dP 

BF  p_  m   wan   g,p 

2    0  dz     <J>    dx 


P  (8-0)  =  P   . 

P         P,i 


dP  (z-X) 

— ^ 0 


(3.1-4) 

(3.1-6) 
(3.1-7) 
(3.1-7) 

(3.1-13) 
(3.1-14) 

C3.1-15) 
(3.2-7) 

(3.2-8) 
(3.2-9) 

(3.2-10) 


b,a    b,v 


(P  ,  "  P   )  - 
P.i    P.e 


(3.2-16) 


IV.   SOLUTION  OF  DIFFERENTIAL  EQUATIONS 

The  model  presented  here  consists  of  solutions  to  the  set  of  equations 

appearing  in  Table  3-1.   Five  particular  cases  of  the  model  are  considered. 

In  the  Cases  1,  2,  3,  and  4  it  is  assumed  that  CL  and  P.  are  connected 

b      b 

by  a  relationship  of  the  form 


where  g  and  g  are  constants  independent  of  variations  in  either  P,  or  x. 

b 

Only  one  gas  component  is  treated. 

The  differences  between  these  four  cases  lie  in  the  values  assigned 
to  the  diffusion  constants.  A  diffusion  constant  equal  to  infinity  implies 
that  diffusion  occurs  so  easily  that  no  partial  pressure  gradient  can  be 
maintained;  a  diffusion  constant  equal  to  zero  implies  no  diffusion;  a 
finite,  non-zero  diffusion  constant  implies  diffusion  as  it  actually  occurs. 
Diffusion  is  assumed  to  occur  only  in  the  gas  phase;  the  diffusion  constant 
of  blood  is  taken  equal  to  zero  in  all  the  models.   Diffusion  can  occur 
from  gas  to  blood  or  from  blood  to  gas,  but  not  within  the  blood  itself. 
The  diffusion  constants  of  air  capillary  gas  and  parabronchial  gas,  n  and 
8  respectively,  are  assumed  different  in  the  derivation  of  Case  4  to  allow 
for  the  possible  effect  of  turbulence  in  the  parabronchial  gas  on  the 
effective  diffusion  constant. 

In  Case  1:  n  =  °°;  8  =  0 

In  Case  2:   n  =  ■»;  8  is  finite,  none-zero 

In  Case  3:   n.  is  finite,  non-zero;  8  =  0 

In  case  4:  n   is  finite,  non-zero;  8  is  finite,  non-zero 
It  can  be  shown  that  the  Case  1  treated  here  is  equivalent  to  the  model 
of  Scheid  and  Piiper  [2]. 
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4.1  Case  1  n  =  ">;  8  =  0 

The  differential  equation  describing  partial  pressure  of  the  gas  com- 
ponent in  the  air  capillary  is  (3.1-13)  subject  to  the  boundary  conditions 
(3.1-14)  and  (3.1-15).  Letting  n  +  -  causes  (3.1-13)  to  take  the  form 


—f  -  0  .  (4.1-1) 

dx 


This  implies 

dP 
dx 


»  constant  .  (4.1-2) 


Equations  (4.1-2)  and  (3.1-15)  imply 

dP 

d^=0   •  W-1"3) 

This  implies  that  P  ,  as  a  function  of  x,  is  a  constant.   Equation  (3.1-14) 

then  implies  that  along  the  entire  length  of  the  air  capillary  the  partial 

pressure  equals  that  of  the  parabronchial  air  the  capillary  opens  into. 

Thus, 

P  =  P    =  P  (4.1-4) 

g    g,P    P  V     ' 

The  differential  equation  describing  partial  pressure  in  the  blood 

capillary  is  (3.1-6)  subject  to  the  boundary  conditions  (3.1-7)  and  to 

(4.1-4).   The  solution  is 

pb "  Pb,v  «•#  +  V1  -  exp(fr»    ■  C4.1-5) 

Equation  (4.1-5)  gives  the  partial  pressure  of  the  gas  component  in  the 
blood  capillary  as  a  function  of  location  along  that  capillary. 

The  differential  equation  describing  partial  pressure  of  the  gas  in 
the  parabronchus  is  (3.2-7)  subject  to  the  boundary  conditions  (3.2-9) 
and  (3.2-10).   Setting  9=0  causes  (3.2-7)  to  take  the  form 


21 


dp      Of 

_E  =  WE£  (P   _  p   }  (4.1-6) 

dz    BF   k  bv    b,a;  l  L   0> 

To  obtain  the  right  side  of  (4.1-6),  x  is  set  equal  to  one  in  (4.1-5).   This 

leads  to 

(Pb,v  "  Pb,a>  +  (Pb,v  "  Pp)(1  "  e*P#»  •      <*•**> 
Equation  (4.1-7)  substituted  into  (4.1-6)  yields 

d^  -  W-  a  -  exp(fr))(pb,v  -  V  •  ^■1"8> 


A  constant,  H,  is  defined  as 


H  =  w6f(l  -  exp(^p))   .  (4.1-9) 


Equation  (4.1-8)  becomes 


dP 

— 2-  +  —  P  =  —  P  (4.1-10) 

dz    BF  p   BF  b,v   '  K     L     V) 

The  solution  of  (4.1-10)  subject  to  (3.2-9)  is 

PP  ■  pP,i  exp(Ii> +  pb,v(1  -  ^O  •        t4-1-11) 

Equation  (4.1-11)  gives  the  partial  pressure  of  the  gas  in  the  parabronchus 
as  a  function  of  location  along  the  parabronchus.  The  partial  pressure 
of  the  gas  component  in  the  arterial  blood  leaving  the  parabronchus  is 
obtained  from  (3.2-16),  recognizing  that  9=0  in  that  expression.  The 
quantity  P    required  in  (3.2-16)  is  found  by  replacing  z  with  X   in 
(4.1-11).   These  substitutions  yield 

Pb,a  *  Yv^  "  f  CI  -  expC^))]  +  F    f-  0.  -   exp(^§))   . 


(4.1-12) 
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The  form  of  (4.1-10)  Is  important.   It  Implies  that  the  manner  in  which 
a  substance  diffuses  from  the  venous  blood  through  the  intervening  blood 
capillary-air  capillary  structure  into  the  parabronchial  gas  is  equivalent 
to  the  manner  in  which  the  substance  would  diffuse  from  a  homogeneous 
compartment  of  venous  blood  across  a  simple  interface  into  the  parabronchial 
gas  if  the  interface  were  characterized  by  a  diffusing  capacity,  H.   Con- 
ceptually, the  blood  capillary-air  capillary  exchange  structure  can  be 
replaced  by  an  equivalent  interface  between  venous  blood  and  parabronchial 
air  compartments.  This  replacement  is  valid  only  insofar  as  the  assumptions 
of  Case  1  are  valid. 

4.2  Case  2  t|  »  ■;  9  is  finite,  non-zero 

In  Case  2,  as  in  Case  1,  n.  =  ™.   Hence  the  same  solutions  are  obtained 
for  P  and  P.  along  the  paired  capillaries.   They  are 


P     =  P.        =  P      , 
g         b.P         P 

and 

Pb   "  Pb,v  -»«tI>   ■ 

Also, 

b,v          h,a            b,v 

Pp(l  -  exp(^f) 


(4.2-1) 
(4.2-2) 


Pp)(l  -  exp(^))  .        (4.2-3) 

The  differential  equation  describing  partial  pressure  of  the  gas 
component  in  the  parabronchus  is  (3.2-7)  subject  to  the  boundary  conditions 
(3.2-9)  and  (3.2-10).   Substituting  (4.2-3)  into  C3.2-7)  yields 

d2pp   BF^p.   .         H_ 

dz2    8<{>  dz  "  6$  p   "  8$  Pb,v  "  (4.2-4) 

where  H  is  defined  by  (4.1-9).   i|i  and  t  are  defined  as 


*  mh     >  and      *"!!••  C4.2-5) 
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Now  the  differential  equation  appears  as 


d  P  dP 

— r2-  -  I  t1  -  i|iP  =  -*P,    •            (4.2-6) 

,2  dz      p      bv 

dz  r 


An  assumed  solution  of  the  homogeneous  equation  is 


P  =  e"   ■  (4.2-7) 

P 


Substitution  into  the  homogeneous  equation  yields 


s2  -  ts  -  i|i  =  0  .  .  (4.2-8) 


which  has  the  roots 


*x   -  J   (t  +^t2  +  4*)   ,  (4.2-9) 


and  s2  =  j   (t  -"\/t2  +  4*)  .  (4.2-10) 

The  solution  to  (4.2-6)  is  of  the  form 

pp  =  v 1  +  v 2  +  pbv  >  (4-2"u) 

be  determined. 
The  boundary  condition  (3.2-9)  applied  to  (4.2-11)  yields 

Pp,i  "  bl  +  b2  +  Pbv   •  (4-2"12> 

The  boundary  condition  (3.2-10)  applied  to  (4.2-11)  yields 

s  A        s9x 
0  ■  bjSj*    +  b2s2e  -    .  (4.2-13) 

Equations  (4.2-12)  and  (4.2-13)  can  be  solved  for  b  and  b„,  yielding 

h'—^r — i^(pP,i-pbv)  -  <4-2-14' 

s2e    -  Sle 
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and 


s,X   p,i    bv 


(4.2-15) 


Substitution  of  (4.2-14)  and  (4.2-15)  into  (4.2-11)  yields  the  solution 

S-A  s.z        S.X  S-Z 


p  =  p,  +  (p  .  -  p.  ) 

p    bv     p,i    bv 


•1" 


(4.2-16) 


The  derivative  with  respect  to  z  is 


dP 

_E 

dz 


(Yi '  V 


s2X  s  z    s  A  s2z 
s1s2(e   e    -  e   e   ) 

s_X      s,>. 
s.e    -  s^ 


(4.2-17) 


The  partial  pressure  of  the  dissolved  gas  component  in  the  arterial  blood 
is  obtained  from  (3.2-16).   The  P    appearing  in  (3.2-16)  is  obtained  by 
letting  z=X  in  (4.2-16).  The  dP  /dz  appearing  in  (3.2-16)  is  obtained  by 
letting  z=0  in  (4.2-17). 


P.    -  P,   +  (P   .  -  Pu  )  & 
b,a    bv     p,i    bv   |  ■ 


5, A      S-jAl 
-  sxe   J 


A)  ■ 


(4.2-18) 


In  the  development  of  Case  2  given  above  it  is  assumed  that  the  partial 
pressure  at  the  entrance  to  the  parabronchus  (z=0)  is  known.   In  an  actual 
experimental  situation,  the  partial  pressure  at  exactly  the  entrance  point 
may  not  be  measurable;  the  partial  pressure  which  might  be  measurable  is 
that  of  the  gas  some  distance  upstream  from  the  entrance.   This  upstream 
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gas  could  be  tracheal  gas  or  air  sac  gas  in  the  real  bird.   Because  of 
diffusion  occurring  between  the  gas  in  the  parabronchial  lumen  and  the 
gas  which  has  not  quite  yet  reached  the  parabronchus ,  the  partial  pressure 
of  a  gas  component  at  the  parabronchial  entrance  will  not  be  the  same  as 
its  partial  pressure  upstream.  This  effect  will  be  greater  at  slower  gas 
flow  rates. 

The  geometrical  relationships  among  preparahronchial,  parabronchial, 
and  postparabronchial  gases  are  taken  to  be  represented  as  shown  in  Figure 
4.  Qualitatively,  the  manner  in  which  the  partial  pressure  of  a  single 
gas  component  (CO.)  changes  is  also  shown.   Although  Figure  4  does  not 
accurately  represent  the  geometrical  relation  between  secondary  bronchi 
and  parabronchi,  it  should  provide  some  insight  into  the  diffusion  processes 
occurring  in  the  neighborhood  of  the  parabronchial  entrance. 

Let  P  ,  P    and  P__  denote  the  partial  pressures  in  the  three 
regions,  respectively.  The  boundary  conditions  of  this  problem  are: 


PT  =  p  .  , 
I   p,l. 

and 

dz  I 

(4.2-19) 

p  =  p 

I   II 

and 

dz  I   dz  II 

(4.2-2Q) 

dz  II 

(4.2-21) 

5  pm  "  °  ■ 

(4.2-22) 

at  z>\    :     -t— 
dz 

The  differential  equation  for  region  I  is  found  by  setting  war)  =  0  in 

(3.2-8),  which  yields 

A  dP 

-2--T^-=0,  (4,2-23) 
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BLOOD  AND   AIR 
CAPILLARIES 


VENOUS  BLOOD 


REGION    I 


REGION     II 


REGION    III 


Figure  4.      A:    Geometrical   relationship  of  pre-parabronchial ,   parabronchial , 
and  post-parabronchial    regions.      B:    Qualitative  nature  of  CO2  partial 
pressure  along   the  parabronchus . 
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where  again 


A  solution, 


(4.2-24) 


Pj  -  ea'  ,  (4.2-25) 

is  assumed.   Substitution  into  (4.2-23)  yields 

s(s-t)  =  0  ,  (4.2-26) 

which  has  the  roots 

s=0  ,   and  s=t  .  (4.2-27) 

The  solution  to  (4.2-23)  is 


(4.2-28) 


Condition  (4.2-19)  allows  b.  to  be  determined.   Then, 
4  ■ 


PT  =  b.e   +  P   .  .  (4.2-29) 

I    3      p,i 

The  general  expression  for  partial  pressure  of  parabronchial  gas  in  region 

II  which  has  already  been  found  is  given  by  (4.2-11): 

pu =  v*1* +  v*2* +  pb,v  •  <4-2-30> 

Next  b  ,  b,,  and  b_,  must  be  found.   From  (4.2-29)  and  (4.2-30)  and  the 
boundary  conditions  (4.2-20)  and  (4.2-21),  the  following  are  determined: 

b3  +  PP,i  "  fal  +  b2  +  Pbv  '  <4-2"31) 

(4.2-32) 

s.X        s.X 
and  0  =  bjSje    +  b2s2e    .  (4.2-33) 

The  three  equations  can  be  solved  for  the  b's: 

s.X 
xs.e   (P   .  -  p.  ) 

bi " v     s2x    P^    s;i i-r  .      ^-^ 

SlS2^e    "  e   '  +  T^s2e    "  sle   ) 
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s.X 
„  ~TSle        (Pp,l  -  Phv> 

b2 ^ — rx — r* rr  •      (4-2-35> 

s1s2(e         -  e       )  +  t(s2e         -  s^e       ) 

s,X  s   X 

s  s,(e         -  e       )(P      .   -  P,    ) 

and  b     =  ^-? Ej^ *SL  (a   j   -,(,<, 

3  Sj_X         s2X  s2X  s  X        •  C4.2-36) 

SlS2^e  "   e        '   +  T(s?e  "  sie        ) 

Substitution  of    (4.2-34),    (4.2-35),    and    (4.2-36)   into    (4.2-29)   and   (4.2-30) 
yields 

tz,  S1X    s2\ 
s.s.e  (e    -  e   ) 

p  =  p    +  (p    -Pi  — 

I    p.i   k  p,i  "W  ■  A    s2X        s  X      s  X   ' 

SlS2^e    "  e   '  +  T^s2e    "  s1e   * 

(4.2-37) 

and  s^  s2X      s2z  s.X 

ts.e   e    -  ts.e   e 


II    bv     p,i    hv       s  X    s.X        s.X      s  X 
s^e    -  e   )  +  T(s2e    -  s^  L   ) 

(4.2-38) 

To  find  the  partial  pressure  in  the  arterial  blood,  substitute  into  3.2-16. 

P„   is  given  by  (4.2-38)  with  z=X;  the  derivative  dP   . /dz  is  taken  equal 
Pie  p,i  n 

to  zero,  since  for  z  much  less  than  zero  the  slope  of  the  parial  pressure 
curve  approaches  zero. 

S1X  S2X  , 

P.    =  P   +  (P   .  -  P,  )  —  1 ; : 2- i . 

b»a    bv     p,L    bv'   •  s  X    s.X         s.X 


"X      2  2        1 

SlS2(e    "  e   ^  +  T(s2e    "  sie   ) 


(4.2-39) 
Equations  (4.2-37)  and  (4.2-38)  give  the  partial  pressure  along  the  parabronchus 
and  (4.2-39)  gives  the  partial  pressure  in  arterial  blood  when  the  boundary 
conditions  specify  the  partial  pressure  of  gas  at  a  location  far  upstream 
(z-"-«>). 
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4.3  Case  3  n  is  finite,  non-zero;  9=0 

In  case  3  diffusion  within  the  air  capillary  is  taken  into  account. 
The  differential  equation  describing  the  partial  pressures  P  and  P  along 
the  paired  capillaries  are  (3.1-6)  and  (3.1-13),  subject  to  the  boundary 
conditions  (3.1-7),  (3.1-14)  and  (3.1-15). 

Figure  5  is  a  simulation  diagram  of  (3.1-6)  and  (3.1-13). 

The  state  variables  are  selected  as 


»1   'l    '    X2  " 

Quantities  a   and  p  are  defined  as 


h 

a  =  — 
n. 


The  state  variable  equation  is 


and 


r,  -  \ 


L'3J 


and   p  =  -jj- 


I- '3 


Taking  the  Laplace  transform  of  (4.3-3)  and  rearranging  yields 


-1 


0 


o 


rv 

y10 

Y2 

" 

y20 

LY3J 

-y30  . 

(4.3-1) 


(4.3-2) 


(4.3-3) 


(4.3-4) 


where  s  is  the  transform  variable;  the  Y's  are  Laplace  transforms  of  the  y's; 
and  the  subscript  0  on  the  y's  signifies  their  initial  values  at  x=0. 
Multiplying  both  sides  of  (4.3-4)  by  the  inverse  of  the  matrix  appearing 
there  yields 


s  (s+p  ) 


(s+p) 
s(s+p) 


-so 
2 


1 


10 


30  J 


(4.3-5) 


30 


31 


where 


-,r"2 

-p   Vp  +  4g 


ands2=^-VpJ-ti2 


(4.3-6) 
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The  technique  of  expansion  in  partial  fractions,  when  applied  to  (4.3-5), 
yields 

r 


|  M2   ,  M3 
s-s    s-s. 


l30 


(4.3-7) 


where 


0   p   -a 


M.  =  — —  i  0   0 
0   p 


(4.3-8) 


s^s^p)   S]+p 


s1(s1+p)   -s^ 

2 
p         s  -a 


(4.3-9) 


and 


s2(s2+p)    s2+p 
s.o 


The  solution  for  the  state  variables  is 


(4.3-10) 


yll  s,x      .,,11' 

=  M  +  Me    +  M  e    [  |  y. 


L  JJ    L 

which  can  be  written  as 


I  '10 

I  '20 

^30 


m, ,        m. 


11        m12        "13 


«,J    |J 


20 


y3i       ^m31       m32       m33j[y30j 


(4.3-11) 


(4.3-12) 


where 


(s.+p)        s..x        (s7+p)        s.x 

e        +  7TTTT  e  •  C4.3-13) 


S.+p  S    X  S.+p  s.x 

+  .    Crr'         +.#*!•  .  (4.3-14) 


■t,  -  3*-*  — T^ r  e8^  +  —7=2 r  e"2*     ,  (4. 3-15) 


(s^+p)        s^        (s,+p)        s.x 


(4.3-16) 


e         +   r.  _.    ■,    e  •  (4.3-17) 


s.x  s.x 

e         +   r/L   S   e  -  (4.3-18) 


S1X                n  s?x 

e   x     +  t—2 r  e  ,  (4.3-19) 


s.x  s.x 

e         +  — r-£ r  e  .  (4.3-20) 


2  2 

s.    -  a         s5x         s.  -  a         s.x 

*a'*fy  +  TfcZ?'         +^?V6  •  (4-3"21) 

Equation    (4.3-12)    is  not   yet   a  usable  solution.      The  boundary  conditions, 
given  by   (3.1-7),    (3.1-14),   and    (3.1-15),    are  constraints  on  different 
quantities   than  those  which  appear  on   the   right    side  of    (4.3-12).      The 
solution  must   be   given  in  terms  of   the  values   of  y?   and  y,  at   x=0  and  y. 
at  x=£.      Substituting  x=4   into    (4.3-12)    and  denoting  by   the  subscript,    «., 
the  values  taken  at  x=£,   yields 


J3 


"yn 

y2i 

= 

y3i 

'm 

"12JI 

m13i 

211 

m22* 

m23£ 

l3U 

m32* 

m33* 

'10 
y20 
y30 


(4.3-22) 


The  known  boundary  values  are  y 


It*  '20' 


and  y™>  t'le  unknown  boundary  values 
are  y  ,  y.,. ,  and  y,^.   Solving  for  the  unknowns  in  terms  of  the  knowns  yields 


-1 

miu 


-1 

_m2U 
L"m3U 


■12* 

(m21£m122.  "  minr22r         V~2U"*13J!.       ""11Z""23JI 
(m3Um12jl  "  mlUm32Z)        (m3Um13£  "  mlUm33l? 


~yu 

y20 

_y30 

(4.3-23) 


[y 


U   "  m12Jly20 


n13*y30J 


(4.3-24) 


It  follows  from  (4.3-24)  that 


'30 


J. 

hia 


1 

"m12£ 

0 

0 

mlH 
0 

-m13J 


'U 


Substituting  (4.3-25)  into  (4.3-12)  yields 


y! 

1 

y2 

miu 

Ly3. 

mll  '"Vlll  -  nllal2t)  (n13"m  "  mllm13)i) 
m21  (m22mll£  "  "21*12^  (m23mlU  "  n21mm) 
m31   (m32mlU  "  m31,,,121)    (m33mlU  "  m31m13£) 


(4.3-25) 


■  20 
r3fl 

(4.3-26) 

d 


Equation    (4.3-26)    is   the  solution  for  y.  ,   y_,    and  y,    (that   is  P    ,  4~  P    , 

12       3  g'  dx  g' 

and  Pfa)  as  a  function  of  x  along  the  paired  capillaries  when  diffusion 
within  the  air  capillary  is  taken  into  account,  and  it  is  given  in  terms  of 
the  known  boundary  conditions.   It  is  equivalent  to  the  following  three 
equations : 


"11  r   ,  (m13mm  -  mllm13£) 


11£ 


"111 


(4.3-27) 
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dpg   m2i       (ln23mll£  "  m21°W 

-r^-^^-P  +    J  11£ Ul     Pu   ,       (4.3-28) 

dx    "ill  P         »Ut         h'v 

m        (m33mHJ  "  m31m13») 

and  P,  =  -S-  P  +   g  11A J1  iJl     v  .       (4.3-29) 

b   "ill  P         "lit         b>v 

The  differential  equation  describing  partial  pressure  in  the  parabronchus 
is  (3.2-8)  subject  to  the  boundary  conditions  (3.2-9)  and  (3.2-10).  The 
specialization  9=0,  which  implies  there  is  only  bulk,  flow  through  the  para- 
bronchus, causes  (3.2-8)  to  assume  the  form 

dP        dP 

—E.  =  3   SiP  (4  3_30x 

dz     BF    dx  ■  «**3  Ju; 

Substituting  (4.3-28)  into  (4.3-30)  and  rearranging  yields 
^£   -wan  rut)     _   .  fW23*      }„ 

The  quantities  H  and  5  are  defined  as 


H   =  wan  P^l  ,  (4.3-32) 


Then  (4.3-31)  assumes  the  form 

dP    u  u 

F+f^f.^,-  (4-3"34) 

Equation  (4.3-34)  has  the  solution 

PP  "  pP,i  <**[=§}  +  pb,v(1-« f1  "  «pH?))  •    (4-3~35) 

Equation  (4.3-35)  gives  the  partial  pressure  of  the  gas  component  in  the 
parabronchus  as  a  function  of  location  along  the  parabronchus  when  there 
is  diffusion  within  the  air  capillaries  and  only  bulk,  flow  through  the 
parabronchus.   The  observation  made  in  connection  with  Case  1  that  the  blood 
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capillary-air  capillary  structure  can  be  replaced  by  an  equivalent  simple 
interface  between  blood  and  gas  compartments  is  precluded  in  Case  3  by  the 
appearance  of  6^0. 

The  partial  pressure  of  the  dissolved  gas  component  in  the  arterial 
blood  is  obtained  from  (3.2-16),  recognizing  that  8=0  and  that  P    is 
obtained  by  letting  z=\  in  (4-3-35).  The  result  is 

+  Pn  i  <Z7>  CI  -  expA)  .  (4.3-36) 

P'1  8Q  BF 


4.4  Case  4:  n  is  finite,  non-zero;  9  is  finite,  non-zero 

In  Case  4,  as  in  Case  3,  n  is  finite  and  non-zero.   Hence  the  same 

dP 
solutions  are  obtained  for  P  ,  P  ,  and   ,SP  along  the  paired  capillaries. 


dx 


They  are 


mll       ^"lli  "  mllm131) 
8   »m  P        ^  h,v  « 

m31       (m33mlli  "  m31m13£) 

Pb  =  ST*-  p„  +        J£  P.    ,        (4.4-2) 

b   "lU  p  mllJt        h'v 

^£=!2UP  .  (°2Um13.  -  mm"23*>  p 
dX    mll£  »  mlU  b'v 

The  differential  equation  describing  partial  pressure  in  the  para- 
bronchus  is  (3.2-8)  subject  to  the  boundary  conditions  (3.2-9)  and  (3.2-10). 
Substituting  (4.4-3)  into  (3.2-8)  and  rearranging  yields 


H  —tL       (™H  ^21A.  .  _  -waa  ^2U    .  mll£m23£1„ 

9*  dz    %     mmJ  Fp    9*  mlu    Lm13£      i^'Vv 


(4.4-4) 
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1  +  dd  ■ 

:  m13Jl  " 

mlUm23l 
m21l 

(4.4-5) 

BF  wan      211 

rd        6<t>        :        *d  "   8*     mm 


Now  the  differential  equation  appears  as 


— B*  -  X.   3-*  -  *,P  ■  -<l>.  (1+6, )P.   .         (4.4-6) 

,2     d  dz     dp     d    d  bv 
dz  r 

By  comparison  of  (4.4-6)  with-  (4.2-6)  it  is  apparent  that  the  solutions  for 

Case  4  are  obtained  from  the  Case  2  solutions  by  making  the  replacements 

T"Td  ' 
•  *♦«  , 

Pbv  *   (1+VPbv   •  ^4.4-7) 

Sl"Sdl     ' 
and  s2*sd2      - 

where  s   and  s,_  are  functions  of  *.  and  t..  Assuming  that  P   .  at  the 
dl      d2  d      d  p,i 

parabronchial  entrance  is  known  and  substituting  into  the  case  2  solutions 
(4.2-16),  (4.2-17),  and  (4.2-18)  yields 

sd2X  sdlz      sdlX  sd2z 

S   £      6       —  S   ti      Q. 

P   =  (1+6, )P,   +  (P   .  -  (1+6.) P.  )   — : — : . 

P       d  bv     p,i    v   d'   bv'  s,„X       s„X 

d2         dl 

Sd2e     "  sdle 

(4.4-8) 

HP  .     ,  Sd2X  SdlZ    SdlX  Sd2Z. 

n  s,.s.,(e    e     -  e    e    ) 

-rf-  =  (P   .  -  (1+6  .)  P.  )  -a^-Si , . (4.4-9) 

dz      P.i      d  bv'  s.-X      s  X      '    K  ' 

„  d2         dl 
s,  e     -  s.,e 
dz         dl 


and 
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Va  "    <1+VPbv  +    <Pp,i  "    ^d^bvM 


*ar 


,   Sd2X    sdl\  "" 


C4.4-10) 


Bd2e     "  3dle 

When  diffusion  occurs  within  the  air  capillaries,  diffusion  and  bulk  flow 
occur  within  the  parabronchus,  and  the  state  of  the  gas  at  z=Q  is  known, 
then  (4.4-8)  and  (4.4-9)  describe  the  partial  pressure  and  its  derivative 
along  the  parabronchus  and  (4.4-10)  describes  the  partial  pressure  of 
arterial  blood.  Assuming  that  P   .  at  a  location  far  upstream  is  known 
and  substituting  into  (4.2-37),  (4.2-38),  and  (4.2-39)  yields 


PI  *  Pp,i  +  (Pp,i-<1+VPbv> 


Td2,  SdlX    Sd2\ 
s,,,e   (e     -  e    ) 


adl"d2 


sJlSj„(e  dl  -e  d2  )  +  t.Cs^e  d2  -sJne  dl  ) 


'dwd2  = 


°dl= 
(4.4-11) 


PTT  =  (l+fi.)P.   +  (P   .-(1+6.)P.  ) 
II       d  bv     p,i     d  hv 


sdlz  sd2X        sd2z  sdlX 
Td3d2e    e     -  Tdsdle    e 

,  SdlX  Sd2\  ^   ,    sd2X 


'dv-  d2" 


Pb,a=    (1+VPbv+(Pp,i-    <1+VPhv>5  I1 


(4.4-12) 


sdlX  sd2X 
V 1    (sd2  -  »dl> 


SH1X    sa?X  s.J      Sj1X 

««<-     "  *  ^  )  +  Vd2-     -  Sdle  ^ 


dl  d2 


.  (4.4-13) 


>J 
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When  diffusion  occurs  within  Che  air  capillaries,  diffusion  and  bulk  flow 
occur  within  the  parabronchus,  and  the  state  of  the  gas  at  z-»— «>  is  known, 
then  (4.4-11)  and  (4.4-12)  describe  the  partial  pressure  along  the  para- 
bronchus and  (4.4-13)  describes  the  partial  pressure  of  arterial  blood. 
In  APPENDIX  B  an  independent  derivation  is  given  for  the  solutions 
for  partial  pressure  along  the  blood  capillaries,  air  capillaries,  and 
parabronchus  when  it  is  assumed  that  diffusion  occurs  within  the  air 
capillaries,  both  diffusion  and  bulk  flow  occur  within  the  parabronchus, 
and  the  state  of  the  gas  at  z-*-°°  is  known. 

4.5  Case  5 

In  Cases  1,  2,  3,  and  4  it  was  assumed  that  the  partial  pressure,  P, 
and  concentration,  C,  of  a  gas  component  in  blood  were  connected  by  a  linear 
relationship  of  the  form 

C  =  k.  +  k2P  , 

If  the  gas  component  of  interest  is  oxygen  the  actual  situation  is  more 
complicated  than  this.   The  oxygen  dissociation  curve,  a  curve  of  concentra- 
tion vs.  partial  pressure,  appears  as  shown  in  Figure  6. 

Due  to  the  effects  named  after  Bohr  and  Haldane,  there  is  an  interaction 
between  oxygen  and  carbon  dioxide  concentrations.   A  change  in  oxygen  partial 
pressure  will  cause  a  shift  of  the  carbon  dioxide  dissociation  curve,  and 
a  change  in  carbon  dioxide  partial  pressure  will  cause  a  distortion  of 
the  oxygen  dissociation  curve.   In  Case  5,  as  in  the  previous  cases,  inter- 
actions between  oxygen  and  carbon  dioxide  are  neglected  and  a  single 
component,  oxygen,  is  treated.   This  component  is  assumed  to  have  the 
dissociation  curve  appearing  in  Figure  6,  and  the  curve  is  not  distorted 
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Figure  6.   Oxygen  dissociation  curve,  linear  approximation  to 
it,  and  2nd  degree  polynomial  approximation  to  it. 
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due  to  the  Bohr  or  Haldane  effects.   The  technique  of  case  5  has  not  been 
used  for  carbon  dioxide  since  this  gas  has  a  dissociation  curve  that  is  very 
linear  over  the  physiological  range. 

In  Case  5  the  same  assumptions  regarding  the  diffusion  constant  are 
made  as  in  Case  1;  that  is  8=*>  and  q=Q.   It  can  be  seen  from  the  computed 
results  of  Figure  9A  that  this  is  a  reasonable  assumption  if  the  ventilation 
rate  is  not  too  low.  The  same  argument  given  in  the  discussion  of  Case  1 
implies  here  also  that  P  equals  a  constant  along  the  length  of  an  air 
capillary. 

In  Cases  1,  2,  3,  and  4  the  solution  was  found  for  partial  pressure, 
P.  ,  along  the  blood  capillary.  By  contrast,  in  Case  5  it  is  more  convenient 
to  solve  for  the  concentration,  C,  .  The  differential  equation  describing 
gas  diffusion  between  an  air  capillary  and  a  blood  capillary  is  (3.1-4)  subject 
to  the  boundary  condition  (3.1-7).   Knowledge  of  the  dissociation  curve,  as 
given  in  Figure  6,  implies  existence  of  a  mathematical  relationship  giving 
P,  as  a  function  of  (1  .   It  is  assumed  here  that  such  a  function  can  be  well 


Pb  =  12Cb  +  "A  +  %      ■  (4-5"1) 

Substitution  of  (4.5-1)  into  (3.1-4),  rearrangement  and  integration,  yields 

dC,  r 

'   h  dx 


+  ^~   =  0  .         (4.5-2) 


«2cb +  "iS +  <vV    J    f 

The  integral  appearing  in  (4.5-2)  is  found  in  reference  [6].   The  lower 


rl  =  ^vy  -  A   •  (4-5"3> 


There  are  two  cases  to  consider. 
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Then 


2  „  -1  2"2Ch  +  ql 
tan   


S 


+  ^  =  0  , 


Jbv 


S-2^+2^tan 


l(2q2Chv  +  qll    ^^ 


(4.5-4) 


(4.5-5) 


1   .    2q2Cb  +  ql  "  ^ 
log 


«T     2q2Ch  +  qx  +  ^ 


+  ^=  0   , 


(4.5-6) 


/-r^l  +  r2exp(-hx/-r1/f))  -  qxQ  -  r2exp(-hx/-r1/f )) 


(4.5-7) 


where 


2"2Cbv  +  ql 


(4.5-8) 


If  a  more  accurate  representation  of  the  dissociation  curve  is  desired, 
in  the  place  of  (4.5-1)  a  third  order  polynomial  can  be  used.   The  integral 
which  then  appears  in  (4.5-2)  can  be  found  in  integral  tables  (reference 
[6]). 
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Equation  (4.5-5)  or  (4.5-7)  gives  the  concentration  of  the  gas  com- 
ponent dissolved  in  blood  if  the  state  of  incoming  venous  blood  is  known  and 
if  the  state  of  the  parabronchial  gas  into  which  the  air  capillary  opens  is 
known.  Due  to  the  form  of  these  expressions,  however,  it  has  not  been 
possible  to  proceed  as  was  done  with-  the  Cases  1,  2,  3,  and  4  and  obtain  an 
analytical  solution  for  the  partial  pressure  along  the  parabronchus.  To 
obtain  these  parabronchial  partial  pressures,  a  discrete  rather  than  con- 
tinuous representation  of  the  parabronchus  has  been  used. 

The  parabronchus  is  considered  to  be  divided  into  a  string  of  cells 
which  are  numbered:   . . .n-l,n,n+l, —   They  all  have  thickness,  Az,  as  shown 
in  Figure  1C. 

It  is  assumed  that  a  value  is  given  for  the  partial  pressure  of  the  gas 
component  in  cell  n,  P  .  Then  C,   is  computable  from  (4.5-5)  or  (4.5-7). 
Setting  x=l  in  the  expression  for  C,   yields  C.   .   (The  additional  subscript 
n  on  C,  indicates  that  the  capillary  is  located  on  cell  n).  From  (3.1-2a), 
the  rate  at  which  the  component  is  transferred  from  a  single  blood  capillary 
into  the  parabronchus  is 

f(Cbv  -  Cban)  '  <4-5"9> 

If  this  gas  enters  cell  n,  it  will  change  the  value  P   which  has  already 
been  assumed  for  that  cell.   So  let  it  enter  cell  n+1.   The  total  gas 
entering  cell  n+1  from  the  capillaries  physically  located  on  cell  n  is 

(4.5-10) 

The  rate  at  which  the  component  flows  directly  from  cell  n  into  cell  n+1  is 

FCpn  ,  (4.5-11) 

where  Cpn  =  BP?n  .  (4.5-12) 
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Under  steady  state  conditions  the  sum  of  (4.5-10)  and  (4.5-11)  equals  the 
rate  at  which  the  component  flows  out  of  cell  n+1.   Thus, 


fw(C   -  C   )Az  +  FC   =  FC    .,  .  (4.5-13) 

bv    ban       pn     p,n+l 


It  follows  that 


fw(C.   -  C,  )Az   +  FC 

P    ,  =  feS tan pn 

p,n+l  BF  ^-0  1,; 

The  process  can  now  be  repeated  to  calculate  P    ,.   In  this  way  the 

partial  pressure  in  all  the  cells  along  the  parabronchus  can  be  calculated, 

given  that  the  partial  pressure  in  the  first  cell  is  P  . . 

The  oxygen  uptake  affected  by  a  single  parabronchus  can  he  expressed 

either  as 

Q(Cba  -  C^)    ,  (4.5-15) 


FB(pPi  -  V  •  (4-5-16> 

Equating  (4.5-15)  and  (4.5-16)  and  solving  for  C,   yields 

-  FR 

C.  .  =  C,   +  —  (P   -  P   )  .  (4.5-17) 

ba    bv    •   pi    pe 
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V.   PARAMETERS  OF  THE  MODEL 

This  section  assigns  values  to  the  various  parameters  of  the  model.  The 
selected  values  have  been  "normalized"  to  fit  a  1.6  kg  bird.  The  description 
is  of  a  female  chicken  (Gallus  Domesticus)  whenever  possible  and  of  a  duck 
(Cairina  Moschata)  when  the  chicken  data  is  unavailable.  Thus  the  avian 
lung  which  is  ultimately  modeled  is  that  of  a  hypothetical,  1.6  kg  female 
bird  resembling  a  chicken  or  duck. 


5. 1  Physical  Parameters  of  the  Lung 

The  parameters  which  define  the  model  parabronchus  are:  I,    X,   a,  4,  w. 
It  has  been  necessary  to  infer  these  parameters  from  descriptions  of  the 
avian  lung  found  in  the  literature  which  have  been  given  in  different  terms. 

The  fractional  distribution  of  lung  volume  for  domestic  fowl  is  given 
by  Duncker  [4]  as  shown  below  in  Figure  7. 
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Figure   7.      Fractional   distribution  of   lung  volume    in  domestic    fowl 
(from  Duncker,   1972) 
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Dunker  also  states  that  " the  air  content  of  the  lungs  is  estimated  by 

taking  the  volume  of  the  air  conducting  parts  plus  50%  of  the  blood  air 
capillary  net  volume  [4]."  Hence,  the  air  capillaries  are  taken  to  compose 
50%  of  the  volume  of  the  parabronchial  mantle. 

Burton  and  Smith  report  that  the  air  volume  of  both  lungs  of  a  1.7  kg 
adult  female  chicken  is  14.6  ml  as  determined  by  buoyancy  on  a  recently 
killed  specimen  [7],  The  adjusted  value  used  here  which  is  appropriate  for 
a  1.6  kg  bird  is  13.74  ml. 

From  the  data  mentioned  above,  the  total  lung  volume  is  determined  to  be 

Vlung, total  =  [(.374)(82.5)  +  (13.9)  +  (. 5) (.566) (82. 5) ]  13"74  ml  =  20-18  ml  ■ 
The  volume  of  all  the  parabronchi  is 

Vrabronchi  "  (-825)  (20.18  ml)  -  16.64  ml   . 
King  and  Cowie  report  that  there  may  be  300-500  parabronchi  in  the  lung 
of  a  chicken  [8].   400  parabronchi  per  lung  are  assumed  here.   The  volume  of 
a  single  parabronchus  is 

V    ,     .    =  16innml  =  .0208  ml   . 
parabronchus      800 

The  length  of  a  parabronchus  in  a  chicken  or  duck  has  not  been  found 
in  the  literature.   But  Dunker  reports  that  the  length  of  paleopulmonic 
parabronchi  of  a  mute  swan  is  20-21  mm  [5].   Based  on  this  fact  a  length  of 
1.7  cm  has  been  somewhat  arbitrarily  selected  for  the  parabronchi  of  the 
model . 

Although,  cut  in  cross-section  perpendicular  to  the  axis  a  parabronchus 
appears  hexagonal,  for  convenience,  the  parabronchi  are  here  assumed  circular 
in  cross-section.  The  diameter  of  a  parabronchus,  d  ,  including  both  the 
mantle  and  lumen  is  found  from 
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So,  the  parabronchial  diameter  is  0.1248  cm.   The  diameter  of  the  parabronchial 
lumen,  which  is  taken  to  include  the  atria,  is  found  from 

d  „l2 

_pj. 


(1.7)  =  (.0208)(.374)  . 

So,  the  parabronchial  lumen  diameter  is  0.763  cm.   The  cross-sectional  area 

2 
of  the  parabronchial  lumen  is  0.C0458  cm  .   The  mantle  thickness  is  0.0242  cm. 

Air  and  blood  capillaries  penetrate  the  entire  thickness  of  the  mantle. 
Also,  speaking  of  the  blood  capillaries  in  the  mantle,  Dunker  states 
"...  the  length  of  the  capillary  stretch  involved  in  exchange  corresponds 
approximately  to  the  thickness  of  the  mantle  of  respiratory  tissue  of  one 
parabronchus  [4]."  Hence,  in  the  model  presented  here  the  length  of  the 
paired  capillaries  is  taken  to  be  0.0242  cm. 

Dunker  reports  the  air  capillary  diameter  in  a  domestic  hen  to  he  0.001 

cm  [4].  The  cross-sectional  area  of  such  an  air  capillary  is  7.854  x  10 

cm  and  the  volume  is  1.901  x  10   cm  .   From  the  total  lung  volume  occupied 

by  air  capillaries,  and  the  volume  of  a  single  air  capillary,  the  total 

number  of  air  capillaries  in  both  lungs  is  found  to  be 

(.5)  (.566)  (16 .64  ml)  8    .„   . 

S =  2.4a  x  10  capillaries 

1.901  x  10 

The  number  of  air  capillaries  associated  with  each  parabronchus  is 
- '- — jJQQ =  3.10  x  10  capillaries  per  parabronchus 


The  linear  density  of  air  capillaries  along  a  parabronchus  is 
•  =  -f-^j ■  1.82  x  10  capillaries  per  cm 
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In  summary,  the  physical  parameters  of  the  model  parabronchus  are  taken 
to  be 

I  =  0.0242  cm   , 

X  -  1.7  cm   . 


*  ■  4.580  x  10"3  cm2 


5.2  Blood  Parameters 

The  0.  capacity  of  avian  blood  used  in  this  model  is  7.2  mmoi'L   , 

which  results  from  a  Hb  density  of  120  gm*L   and  an  Hb  capacity  for  0_ 

of  0.060  mmol-gm  .   This  value  falls  within  the  range  reported  by  Scheipers 

t  «1.  [9]  of  6.7+1.3  mmol'L   for  duck  blood.   Scheipers  et  al.  [9]  also 

display  the  in  vivo  0.  dissociation  curve  for  duck  blood.   For  Cases  1,  2,  3, 

and  4  of  the  model  it  is  necessary  to  approximate  the  sigmoid  curve  by  a 

straight  line,  the  slope  of  which  is  the  value  of  Sn9-  This  is  done  by 

assuming  the  straight  line  to  pass  through  the  82  torr  point  representing 

arterial  blood  and  the  39  torr  point  representing  venous  blood.  The 

resulting  slope,  or  solubility  coefficient  is 

_  (0.91)  -  O.43)7.2mmol-L"1  -1  ro   -1 

S (82  -  39)  torr °-0804  mmol*L  '  torr 

In  Case  5  it  is  required  that  the  0.  dissociation  curve  be  approximated, 

not  by  a  straight  line,  but  by  a  polynomial  of  second  degree.   This 

polynomial,  giving  partial  pressure  as  a  function  of  concentration,  will 

have  the  general  form 

„2 
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The  constants  qQ,  q  ,  q2,  which  are  to  be  chosen  so  that  over  the  physiologi- 
cal range  the  polynomial  accurately  represents  the  dissociation  curve,  have 
been  found  by  picking  three  points  on  the  actual  dissociation  curve, 

1)  Px  =  81  torr    ,   C.  -  6.55  mmol-L-1 

2)  P2  =  55  C2  =  5.256 

3)  P3  =  42  C2  =  3.60  . 

Substituting  them  into  the  general  polynomial  expression,  obtaining  three 
unknowns,  and  solving  the  equations  for  q1 ,  q,,  q  yields 
q„  =  92.29  torr   . 

q,  =  -28.92  torr-L-mmol"1  , 
and  q2 

At  the  three  chosen  points  the  polynomial  approximation  will  be  exact; 
elsewhere  in  the  physiological  range  it  will  be  close  to  the  real  dissocia- 
tion curve. 

The  effective  in  vivo  dissociation  curve,  a  linear  approximation  to  it, 
and  a  polynomial  approximation  to  it  are  illustrated  in  Figure  6.  The 
linear  approximation  is  a  crude  fit  to  the  real  dissociation  curve  over 
the  physiological  range.   Between  about  42  and  100  torr  the  polynomial 
provides  a  quite  accurate  approximation  to  the  real  dissociation  curve. 
Above  100  torr  the  polynomial  approximation  tends  toward  oxygen  concentra- 
tions which  are  slightly  too  high;  below  42  torr  the  approximation  and  the 
real  curve  diverge  widely. 

The  effective  in  vivo  carbon  dioxide  dissociation  curve  for  duck 
blood  determined  by  Scheipers  et  al.  [9]  appears  to  be  linear  over  the 
physiological  range.  They  give  the  slope  of  that  curve  to  be  g  ,  =  0.27 
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mmol-L   -torr   .   It  has  not  been  attempted  to  find  a  polynomial  approximation 
to  the  CO.  dissociation  curve  for  Case  5  of  the  model. 

5.3  Diffusing  Capacity  of  the  Lung 

By  a  carbon  monoxide  diffusing  capacity  test,  Piiper,  Pfeifer,  and 
Scheid  [10]  determined  the  0.  diffusing  capacity  of  the  lungs  of  1.4  kg 
chickens  to  be  DQ2  =1.0  ml-min"  -torr"  .   The  value  to  be  taken  for  the 
model  of  the  1.6  kg  bird  is  1.14  ml-min"  -torr   .   This  is  in  agreement 
with  the  value  which  Scheid  and  Piiper  [2]  determined  by  use  of  a  model 
similar  to  the  Case  1  developed  here.  From  total  lung  diffusing  capacity, 
total  parabronchial  mantle  air  volume,  the  air  volume  of  a  single  capillary, 
and  the  length  of  an  air  capillary,  it  can  be  determined  that  the  value  of 
h,  diffusing  capacity  per  length  along  a  single  air  capillary,  is 
h=l. 901x10   ml-min  -torr   -cm  .   If  mass  is  expressed  in  moles  this 
becomes  h=0. 849x10   mmol-min  -torr   -cm"  .  The  value  used  here  for 
C02  diffusing  capacity  per  length  along  a  single  air  capillary  is 
1.698  x  10   mmol-min  -torr   'cm  ,  which  is  20  times  larger  than  the 
0„  diffusing  capacity. 

5.4  Gas  Parameters 

The  equation  describing  the  diffusion  of  a  single  gas  component  in  air 

written  in  terms  of  concentration  is 

dM     „  dc 
-r-  =  -aD  -j-     . 
dt       dx 

Written  in  terms  of  partial  pressure  it  is 

dM      dP 

Zt"  =  -an  dx"  ' 

where  n  =  BD  . 
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From  the  Handbook  of  Chemistry  and  Physics  (reference  [11])  are  found 
the  values 


and 


B  is  the  solubility  of  a  gas  component  in  air. 
1     mi 


760  mi-torr 


So, 


Tjnol 


m.in*cm*  torr 


=  4.899  x  10"4  -  mmo1 


C0«  min*cm«torr 

These  are  the  values  assigned  to  the  diffusion  constants  in  the  air  spaces 
of  both  the  parabronchial  lumen  and  air  capillaries. 

5.5  Ventilation  and  Cardiac  Output 

It  is  reported  by  King  and  Payne  [12]  that  the  minute  ventilation  of  a 
resting,  unanesthetized  3.4  kg  female  chicken  is  766  ml/min.  Based  on  this 
data,  the  basal,  resting  value  to  be  used  for  the  1.6  kg  bird  of  this  model 
is  360  ml/min. 

Vogel  and  Sturkie  [13]  give  values  of  cardiac  output  for  female  chickens. 
Based  on  their  results,  the  value  to  be  used  here  for  basal,  resting  cardiac 
output  of  a  1.6  kg  bird  is  234  ml/min. 
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VI.   RESULTS  AND  DISCUSSION 

In  Chapter  III  differential  equations  have  been  derived  which  describe 
the  partial  pressures  of  a  diffusing  gas  component  in  the  air  capillaries, 
blood  capillaries,  and  lumen  of  a  model  of  the  avian  parabronchus .   In 
Chapter  IV  solutions  for  these  equations  have  been  found.   Evaluations  of 
those  solutions  which  will  describe  the  performance  of  the  model  have  been 
made,  and  they  are  presented  graphically  in  this  Chapter. 

Figures  8A  and  8B  illustrate  the  partial  pressures  of  oxygen  along  a 
blood  capillary-air  capillary  pair  which  have  been  computed  according  to 
(4.3-27)  and  (4.3-29).  A  linear  dissociation  curve  has  been  assumed  and 
a  finite  diffusion  coefficient  has  been  used  for  the  gaseous  medium  in  the 
air  capillary,  these  parameters  taking  the  values  which  were  arrived  at  in 
Chapter  V.   The  air  capillary  has  been  assumed  to  open  into  parabronchial 
gas  with  P  _  =  160  torr,  this  value  giving  a  relatively  great  partial 
pressure  difference  between  the  blood  and  gas  phases.   Two  cases  have  been 
considered  for  the  state  of  the  venous  blood:   P  „  =  20  torr  and  P._  =  40 
mmHg. 

For  the  dotted  curves  of  Figures  8A  and  8B  the  cardiac  output  has  been 
taken  to  be  Q  «  234  ml/min,  corresponding  to  the  basal,  resting  rate.   For 
both  venous  blood  P  «  values  the  end-capillary  blood  reaches  an  oxygen 
partial  pressure  of  about  150  torr,  coming  to  within  10  torr  0f  complete 
equilibration  with  air  capillary  gas.   In  both  cases  the  P   in  the  air 
capillary  varies  only  0.3  torr  along  its  entire  length. 

For  the  solid  curves  of  Figures  8A  and  8B  the  cardiac  output  has  heen 
increased  to  2340  ml/min.   For  both  venous  blood  P   values  the  blood  P 
rises  only  about  30  torr  along  the  length  of  the  capillary,  increasing  its 
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Figure  8.   Partial  pressures  along  a  blood  capillary  and  an  air  capillary  when 
cardiac  output  is  234  ml/min  (dotted  lines)  and  when  cardiac  output  is  2340  ml/min 


Pvo2*40'  pVC02*45,and  PVC02=55  torr'   The  luminal  Eas  at  the  open  end  of  the  air 
capillary  is  described  by  PQ2«160  and  PCo,-0  torr.   A:   PQ,  along  the  blood 
capillary.   B:   PQ2  along  the  air  capillary.   C:   PCQ2  along  the  blood  capillary. 
D:   Prj02  alonS  the   air  capillary. 
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partial  pressure  only  about  25%  of  the  way  to  full  equilibration  with  air 
capillary  gas.   Since  there  is  a  greater  difference  between  air  capillary 
partial  pressure  and  blood  capillary  partial  pressure  along  the  entire 
length  of  the  capillary  pair,  a  greater  amount  of  oxygen  diffuses  into  the 
blood  from  the  air  than  in  the  case  of  Q  =  234  ml/min.   The  increased  rate 
of  oxygen  diffusing  into  the  air  capillary  then  causes  an  increased  oxygen 
partial  pressure  gradient  along  the  length  of  the  air  capillary,  causing  a 
variation  in  oxygen  partial  pressure  of  about  0.6  torr  along  its  entire 
length.   Thus,  increasing  the  rate  of  blood  flow  through  the  capillary  by  a 
factor  of  ten  increases  both  the  oxygen  partial  pressure  gradient  in  the 
air  capillary  by  a  factor  of  approximately  two  and  the  rate  of  transfer  of 
oxygen  from  the  air  capillary  to  the  blood  capillary  by  a  factor  of  about 
two. 

Figures  8C  and  8D  illustrate  the  partial  pressures  of  carbon  dioxide 
along  the  blood  capillary-air  capillary  pair.  Again,  a  linear  dissociation 
curve  and  a  finite  diffusion  coefficient  have  been  used.  The  air  capillary 
has  been  assumed  to  open  into  parabronchial  gas  with  T?rri7   ■  0.0  torr;  the 


(dotted  curves)  and  Q  =  2340  ml/min  (solid  curves).  Observations  similar 
to  those  made  above  for  oxygen  can  also  be  made  for  carbon  dioxide.   Equili- 
bration of  hlood  with  the  gas  phase  occurs  more  readily  at  the  lower  cardiac 
output;  the  variation  in  carbon  dioxide  partial  pressure  along  the  length 
of  the  air  capillary  remains  below  1.0  torr;  and  increasing  the  cardiac 
output  by  a  factor  of  ten  increases  the  air  capillary  partial  gradient  and 
rate  of  carbon  dioxide  transfer  from  blood  to  gas,  both  by  a  factor  of 
approximately  two. 
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In  a  brief  analysis  of  diffusion  in  the  air  capillaries  of  a  crow's 
lung,  Hazlehoff  [14]  calculated  that  if  the  oxygen  uptake  of  the  bird  is 
2200  ml/min,  the  oxygen  partial  pressure  difference  between  the  ends  of  the 
capillary  is  about  19.0  torr.  The  oxygen  uptake  of  his  example  is  about  84 
times  that  of  the  resting  bird  used  in  this  research  and  his  calculated 
partial  pressure  difference  divided  by  84  is  .23  torr.   This  is  a  value 
very  close  to  that  arrived  at  here  and  illustrated  in  Figure  8B. 

The  very  small  variation  in  oxygen  and  carbon  dioxide  partial  pressures 
along  the  air  capillary  implies  that  with  only  slight  error  the  partial 
pressure  of  these  gases  can  be  regarded  as  essentially  constant  along  the 
air  capillary. 

Figure  9A  illustrates  the  variation  of  oxygen  partial  pressure  along 
the  parabronchus.   Initially,  the  state  of  the  gas  is  described  by  P  ,  = 
160  torr.  As  the  gas  passes  through  the  parabronchus  oxygen  diffuses 
through  the  blood-air  capillary  network  into  the  blood  causing  the  para- 
bronchial  partial  pressure  to  decrease.   Finally,  the  gas  leaves  the  para- 
bronchus with  a  lower  oxygen  partial  pressure  than  it  entered  with.   At 
fixed  cardiac  output  of  234  ml/min  ventilation  rates  from  18  to  720  ml/min 
have  been  assumed  and  the  state  of  venous  blood  is  described  by  P  ,  =  40 
torr. 

Assuming  a  given  average  ventilation  for  the  lungs,  at  any  instant  the 
individual  members  of  the  population  of  parabronchi  in  the  lungs  of  a  real 
bird  will  experience  a  spread  of  ventilatory  rates  both  above  and  below  the 
average  rate.  Also,  the  ventilation  of  the  parabronchi  will  vary  second  by 
second  through  the  different  phases  of  the  respiratory  cycle.  For  these 
reasons  different  gas  flow  rates  through  the  parabronchi  are  of  interest 
even  at  a  constant,  average  whole- lung  ventilation  rate. 
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Figure  9.   Partial  pressures  along  the  parabronchus  computed  under  the  assumption 
of  (1):  equal  diffusion  resistance  in  the  air  capillaries  and  parabronchial  lumen 
(dotted  curves)and,  (2):  constant  partial  pressure  along  the  air  capillaries  and 
bulk  flow  only  in  the  parabronchial  lumen  (solid  curves).   The  ventilation  varies 
from  18  to  720  ml/min.   The  state  of  inspired  gas  is  Pjn2=160  ami  PiC02=0  torr» 
the  state  of  venous  blood  is  Pyo2=4°  and  PVC02=45  torr.   A:   PQ2  along  the 
parabronchus.   B:   Pc02  alon£  tne  parabronchus. 
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For  the  solid  curves  of  Figure  9A  it  has  been  assumed  that  the  gas 
moves  by  bulk  flow  without  any  diffusion  either  forward  or  backward  in  the 
axial  direction.  This  is  the  assumption  of  Case  1,  and  these  curves  have 
been  calculated  according  to  (4.3-11).   It  can  be  seen  that  at  the  slower 
rates,  better  equilibration  of  parabronchial  gas  with  the  venous  blood  is 
achieved.   Another  obvious  characteristic  of  these  curves  is  that  they 
indicate  no  change  in  the  partial  pressure  of  the  gas  until  it  physically 
enters  the  parabronchus. 

For  the  dotted  curves  of  Figure  9A  movement  of  gas  through  the  para- 
bronchus both  by  bulk  flow  and  by  diffusion  in  the  axial  direction  have 
been  accounted  for.   These  results  have  been  calculated  according  to  (4.4-11) 
and  (4.4-12).  The  most  remarkable  characteristic  of  these  curves  is  that 
they  indicate  a  change  in  the  partial  pressure  of  the  gas  before  it  physically 
enters  the  parabronchus.  This  is  due  to  the  gas  which  has  not  yet  reached 
the  parabronchus  diffusing  toward  the  region  of  lowered  oxygen  partial 
pressure  within  the  parabronchus. 

At  high  ventilation  rates  bulk  flow  is  relatively  much  more  important 
than  diffusion,  so  diffusion  can  be  neglected  with  only  a  slight  error 
being  introduced.   At  lower  ventilation  rates  diffusion  becomes  more  impor- 
tant, and  the  partial  pressure  variation  along  the  parabronchus  is  not 
well  approximated  by  the  solid  curves  of  Figure  9A  in  which  diffusion  is 
neglected.   If,  however,  one  is  interested  not  in  the  partial  pressure 
profile  along  the  parabronchus,  but  only  in  the  overall  rate  of  oxygen 
exchange  taking  place  in  the  parabronchus,  it  appears  not  to  matter  very 
much  whether  diffusion  is  accounted  for  or  not.   This  conclusion  follows 
from  a  comparison,  at  a  given  ventilation  rate,  of  the  end-parabronchial 
partial  pressure  values  in  Figures  9A  and  9B  which  reveals  only  a  small 
difference. 
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Figure  9B  illustrates  the  variation  of  carbon  dioxide  partial  pressure 
along  the  parabronchus.   The  gas  approaching  the  parabronchus  is  assumed 
initially  to  have  P    =0.0  torr;  the  venous  blood  is  assumed  to  have 
P„„2  =  ^  torr;  the  cardiac  output  is  assumed  fixed  at  234  ml/min;  and  a 
range  of  ventilation  rates  from  18  to  720  ml/min  is  considered.  For  the 
solid  curves  diffusion  is  neglected;  for  the  dotted  curves  diffusion  is 
taken  into  account.  The  same  conclusions  which  were  arrived  at  in  the 
consideration  of  oxygen  also  apply  to  these  carbon  dioxide  curves. 

The  attempt  is  made  in  the  discussion  of  Figures  10A  and  10B  to  give  a 
more  complete  picture  of  the  whole  parabronchus  as  a  gas  exchange  device  by 
combining  the  results  describing  partial  pressure  variation  along  the  blood 
capillaries  with  those  describing  the  variation  along  the  parabronchial 
lumen.  Figure  10A  treats  oxygen;  Figure  10B  treats  carbon  dioixide.  The 
series  of  curves  on  the  left  describes  partial  pressure  variation  in  blood 
capillaries  which  are  located  at  different  positions  along  the  parabronchus. 
In  each  capillary  the  partial  pressure  changes  from  the  venous  value  to  the 
end-capillary  value,  the  rate  and  extent  of  that  rise  depending  on  the 
partial  pressure  of  the  gas  in  the  adjacent  air  capillary.  The  partial 
pressure  in  the  air  capillary  depends  upon  the  partial  pressure  in  the 
parabronchial  gas  into  which,  the  air  capillary  opens,  which  in  turn  depends 
upon  its  location  along  the  parabronchus.   The  variation  in  partial  pressure 
of  the  gas  along  the  parabronchus  is  shown  by  the  curve  on  the  right.   The 
physical  location  of  a  blood  capillary  with  respect  to  position  along  the 
parabronchus  is  indicated  by  the  dotted  lines.   The  vertical  segment  of  the 
dotted  line  indicates  the  extent  to  which  the  capillary  blood  falls  short 
of  full  equilibration  with  the  gas  in  the  adjacent  air  capillary.   The  gas 
leaving  the  parabronchus  is  at  the  partial  pressure  indicated  by  the  right 
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Figure  10.   A:   O,  partial  pressure  along  several  parabronchial  blood  capillaries 
(left)  and  along  the  parabronchial  lumen  (right)    The  capillary  location  along 
the  parabronchus  is  indicated  by  horizontal  dotted  lines.   The  extent  to  which 
blood  falls  short  of  full  equilibration  with  air  is  indicated  by  vertical  dotted 
lines.   B:   C02  partial  pressure  along  several  parabronchial  blood  capillaries 
and  along  the  parabronchial  lumen. 
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end  point  of  the  parabronchial  gas  curve.   The  arterial  blood  is  formed 
from  the  blood  emptying  out  of  all  the  blood  capillaries  along  the  para- 
bronchus,  so  the  partial  pressure  of  the  arterial  blood  is  some  average  of 
the  values  of  partial  pressure  of  the  right  end  points  of  the  blood  capillary 
curves . 

The  results  illustrated  in  Figures  8  through  10  have  all  been  based 
upon  the  assumption  of  a  linear  dissociation  curve.  This  is  reasonable 
when  considering  carbon  dioxide,  but  not  when  considering  oxygen.  A  con- 
sequence of  the  linear  dissociation  curve  used  in  the  model  so  far  is  that 
when  the  partial  pressure  of  oxygen  in  blood  climbs  above  about  100  torr  , 
the  concentration  of  oxygen  in  blood  rises  above  the  real  physiological 
limits.   Calculations  based  upon  the  model  with  a  linear  dissociation  curve 
can  exaggerate  the  rate  at  which  oxygen  is  transferred  from  the  gas  phase 
to  the  blood  phase.  The  oxygen  concentration  in  the  arterial  blood  can  be 
too  high  and  the  partial  pressure  of  oxygen  in  the  end-parabronchial  gas 
can  be  too  low. 

To  obtain  more  accurate  results  for  oxygen  exchange,  Case  5  of  the 
model  has  been  developed  in  which  the  non-linearity  of  the  oxygen  dissocia- 
tion curve  has  been  taken  into  account.   Based  on  Case  5,  Figure  11  illus- 
trates the  variation  of  oxygen  partial  pressure  in  the  blood  capillaries  at 
different  locations  along  the  parabronchus  and  the  variation  of  oxygen 
partial  pressure  in  the  parabronchial  gas.   The  ventilation  rate  has  been 
taken  to  be  360  ml/min;  the  cardiac  output,  234  ml/min;  the  partial  pressure 
of  venous  blood,  43  torr  .   The  reason  for  increasing  the  partial  pressure 
of  venous  blood  in  Case  5  is  to  keep  on  the  right  limb  of  the  approximate 
dissociation  curve  illustrated  in  Figure  6.   Figure  12,  also  calculated 


60 


(Haoi)  aanssaad  ivuavd  °o 


61 


U     O 

■z.   en 

<     CO 


—  CL 


(yyoi)  Banssaad  ivuavd  Go 


62 


from  Case  5,  is  similar  to  Figure  11,  but  the  ventilation  rate  and  the 
cardiac  output  are  1800  ml/min  and  1170  ml/min,  respectively.   Figures 
13A  and  13B  illustrate  the  oxygen  concentration  in  the  same  blood  capil- 
laries under  the  same  conditions  as  in  Figures  11  and  12,  respectively. 

Comparison  of  Figures  10A  and  11  reveals  some  interesting  differences. 
In  Figure  10A  the  curves  describing  the  rise  of  oxygen  partial  pressure  in 
the  blood  capillaries  have  the  character  of  a  decaying  exponential.   In 
Figure  11  the  blood  capillary  curves  rise  faster  initially  to  a  knee,  then 
form  a  plateau  where  they  are  almost  level.  Also,  in  Figure  11  the  capillary 
blood  appears  to  equilibrate  more  completely  with  the  gas  phase.  The  curve 
in  Figure  11  describing  oxygen  partial  pressure  of  parabronchial  gas  does 
not  fall  to  an  end-parabronchial  value  as  low  as  that  of  Figure  10A.  This 
is  a  direct  consequence  of  the  blood  with  a  more  realistic,  non-linear 
dissociation  curve  becoming  saturated.  Thus,  less  oxygen  is  extracted  from 
the  parabronchial  gas. 

When  ventilation  rate  and  cardiac  output  are  increased  five-fold,  as 
shown  in  Figure  12,  the  capillary  blood  is  quite  far  from  full  equilibra- 
tion with  the  gas  phase  by  the  end-capillary  point  and  the  oxygen  partial 
pressure  of  the  gas  phase  does  not  drop  so  far  by  the  end-parabronchial 
point . 

From  Figure  13  it  can  be  seen  that  at  the  basal,  resting  ventilation 
rate  and  cardiac  output  the  oxygen  concentration  of  capillary  blood  rises 
quickly  to  a  plateau  and  increases  only  slightly  thereafter.   When  the 
flow  rates  are  increased  by  a  factor  of  five  the  oxygen  concentration  in 
the  blood  capillaries  rises  steadily  along  the  entire  capillary  length 
and  reaches  the  end-capillary  point  in  a  nearly,  hut  not  completely, 
saturated  condition. 
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Figure  14  shows  the  oxygen  partial  pressure,  calculated  from  Case  5, 
of  the  end-parabronchial  gas  under  different  conditions  of  ventilation  and 
cardiac  output.   The  inspired  gas  is  assumed  to  have  P  =160  torr ,  and  the 
venous  blood  is  assumed  to  have  P 02=43  torr .  The  trends  are  what  one  would 
expect  intuitively.   Increasing  ventilation  rates  gives  the  gas  less  chance 
to  equilibrate  with  blood,  so  the  end-parabronchial  oxygen  partial  pressure 
increases.  Decreasing  cardiac  output  causes  the  gas  to  be  exposed  to  more 
fully  saturated  blood  on  the  average,  again  causing  end-parabronchial 
oxygen  partial  pressure  to  increase. 

Recall  that  blood  exiting  from  all  the  different  capillaries  in  the 
lung  at  many  different  states  of  saturation  mixes  together  to  form  the 
arterial  blood.  Figures  15  and  16  illustrate  the  arterial  oxygen  partial 
pressure  and  concentration,  respectively,  which  have  been  calculated  accord- 
ing to  Case  5,  under  different  conditions  of  ventilation  and  cardiac  output. 
Increasing  cardiac  output  decreases  saturation  of  arterial  blood,  and 
increasing  ventilation  increases  the  saturation. 

Increasing  cardiac  output  results  in  an  increased  quantity  of 
oxygenated  blood  carried  from  the  lung  to  the  tissues;  however,  increasing 
cardiac  output  causes  that  blood  which  is  carried  to  contain  less  oxygen. 
So,  what  is  the  net  effect?  What  happens  to  oxygen  uptake? 

An  elementary  argument  forces  one  to  conclude  that,  if  the  partial 
pressures  of  the  venous  blood  and  inspired  air  entering  the  lung  are  fixed 
at  constant,  unvarying  values,  then  increasing  the  rate  of  flow  of  either 
blood  or  gas  must  increase  the  rate  of  oxygen  uptake  by  the  blood.   In- 
creasing the  rate  of  flow  of  blood  or  gas  through  the  lung  will  decrease 
the  equilibration  between  blood  and  gas.   This  decreased  equilibration  will 
result  in  greater  partial  pressure  differences  across  the  exchange  surfaces 
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of  the  lung,  which  will  in  turn  cause  a  greater  rate  of  oxygen  transfer 
across  those  exchange  surfaces.  That  is,  greater  oxygen  uptake. 

The  model,  Case  5  with  the  approximate,  non-linear  dissociation  curve, 
has  been  used  to  compute  oxygen  uptake  at  various  values  of  ventilation 
rate  and  cardiac  output.  The  range  is  from  basal,  resting  rates  up  to  eight 
times  normal  for  both  blood  and  gas.   The  inspired  gas  has  P  .=160  torr  and 
the  venous  blood  has  P02=43  torr  .  The  results  are  shown  in  Figure  17 . 

From  Figure  17  it  is  seen  that  the  greatest  increase  in  oxygen  uptake 
is  achieved  by  increasing  both  cardiac  output  and  ventilation  as  much  as 
possible.  An  interesting  result  displayed  by  Figure  17  concerns  small 
increases  in  oxygen  uptake.   Starting  from  basal  flow  rates,  increasing 
the  ventilation  while  keeping  cardiac  output  constant  causes  a  negligible 
increase  in  oxygen  uptake.  However,  increasing  the  cardiac  output  while 
keeping  ventilation  constant  causes  a  very  significant  increase  in  oxygen 
uptake.  This  result  is  explainable.  At  the  basal  flow  rates  blood  that 
leaves  the  lung  is  almost  fully  oxygenated.   Increasing  the  ventilation  will 
not  significantly  increase  the  level  of  oxygenation  of  the  blood,  so  the 
oxygen  uptake  cannot  increase  significantly.   But  increasing  the  cardiac 
output  will  cause  an  increased  amoung  of  blood  that  is  fully  oxygenated, 
or  nearly  so,  to  be  carried  from  the  lung  to  the  tissues. 

By  drawing  a  horizontal  line  across  Figure  17  at  a  particular  level  of 
oxygen  uptake  one  can  see  that  there  is  an  infinite  number  of  combinations 
of  gas  and  blood  flow  rates  that  will  satisfy  that  level  of  oxygen  uptake. 
This  model  does  not  select  any  one  of  those  as  more  preferable  than  the 
others.   It  would  be  interesting  to  determine  experimentally  which  points 
on  such  a  graph  an  exercising  bird  selects. 
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The  maximum  increase  in  oxygen  uptake  shown  in  Figure  17  is  about 
4-1/2  times  above  the  basal  level.  This  occurs  when  blood  and  gas  flow 
rates  are  increased  eight  times  above  their  basal  levels.   This  may  seem 
low.   The  model,  however,  assumes  a  constant  level  of  venous  blood  oxygen 
partial  pressure  regardless  of  cardiac  output  and  ventilation  rate,  which 
is  a  condition  that  probably  does  not  hold  in  a  real  bird.  But  the  venous 
oxygen  partial  pressure  levels  are  determined  by  the  metabolic  rates  of  the 
tissues  and  the  model  constructed  here  is  a  model  of  the  lung  only.   If  the 
venous  blood  were  to  take  a  lower  value  for  oxygen  partial  pressure,  which 
the  Case  5  model  in  its  present  form  cannot  accommodate,  the  oxygen  uptake 
might  be  boosted  toward  a  level  that  better  simulates  exercise  conditions. 

As  ventilation  rate  and  cardiac  output  increase,  the  average  oxygen 
partial  pressure  of  the  gas  exposed  to  the  exchange  surfaces  of  the  lung 
approaches  that  of  fresh,  inspired  air,  and  the  average  partial  pressure  of 
the  blood  exposed  to  those  surfaces  approaches  that  of  venous  blood.   So  an 
upper  limit  on  oxygen  uptake  can  be  calculated  by  assuming  the  exchange 
surfaces  to  be  completely  exposed  to  venous  blood  on  one  side  and  fresh  air 
on  the  other.  The  oxygen  uptake  under  these  hypothetical  conditions  is 
calculated  by  multiplying  the  total  diffusing  capacity  of  the  lungs  times 
the  difference  in  partial  pressure  between  the  blood  and  gas  phases. 
Results  are  shown  in  Figure  18  for  a  range  of  oxygen  partial  pressures  in 
the  blood  and  an  oxygen  partial  pressure  of  160  mmHg  in  the  gas.   The 
diffusing  capacity  has  been  taken  to  be  .0509  mmol-min   -torr"1. 

In  considering  the  oxygen  uptake  increases  shown  in  Figure  17  it 
should  also  be  remembered  that  the  lung  parameters  used  in  this  model  are 
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Figure  18.   Upper  limit  on  O,  uptake  by  the  lungs  as  a  function  of  venous  02 
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those  of  a  female  chicken,  a  bird  which  may  not  be  able  to  achieve  a  high 
oxygen  uptake.   In  particular,  the  air  capillaries  of  other  birds  have  a 
smaller  diameter  than  .001  cm,  reaching  as  low  a  dimension  as  .0003  cm  in 
some  species  (Dunker  [4]).   Smaller  air  capillaries  increase  the  diffusing 
capacity  of  the  lung,  which  in  turn  increase  the  maximum  oxygen  uptake. 

Kawashiro  and  Scheid  [15]  report  that  in  resting,  unanesthetized  hens 
the  oxygen  partial  pressure  of  arterial  blood  is  82  torr.   The  results 
obtained  here  from  the  model  indicate  a  much  higher  arterial  oxygen  partial 
pressure  in  resting  birds.   Since  under  resting  conditions  maximum  oxygena- 
tion of  blood  should  be  obtained  in  the  real  bird  as  well  as  in  the  model, 
this  discrepancy  presents  a  challenge  to  the  assumptions  upon  which  the  model 
is  based.  Two  effects  which  the  model  does  not  take  into  account  are  the 
shunting  of  blood  past  the  lung  without  an  opportunity  to  undergo  gas 
exchange  and  the  presence  of  ventilation-perfusion  imbalances,  whatever  may 
be  the  form  they  assume  in  an  avian  lung. 

In  a  real  bird  some  blood  may  be  shunted  past  the  lung  without  being 
exposed  to  the  gas  exchange  surfaces.  Under  basal  conditions  that  blood 
which  is  not  shunted  will  undergo  gas  exchange  leaving  the  parabronchi 
at  P  _  =  140  torr  and  C  ,  =  8.25  mmol/L.   Its  state  is  described  by  a  point 
far  up  the  dissociation  curve  of  Figure  6.  Mixing  with  a  small  amount 
of  shunted  blood  will  reduce  the  oxygen  concentration  slightly,  but  because 
of  the  slope  of  the  dissociation  curve  the  partial  pressure  will  be  greatly 
reduced,  perhaps  to  82  torr   Ventilation-perfusion  imbalances  should  have 
a  similar  effect. 
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VII.   RECOMMENDATIONS  FOR  FUTURE  WORK 

A  logical  project  with  which  to  follow  up  this  thesis  would  be  an 
experiment  having  the  dual  objectives  both  of  testing  the  accuracy  of  the 
model  and,  if  it  proves  sufficiently  accurate,  of  determining  possible  loca- 
tions of  the  CO.  receptors  in  the  avian  lung. 

To  test  the  model  it  would  be  necessary  to  unidirectionally  ventilate 
the  lung,  as  discussed  in  reference  £3],  of  a  species  having  negligible 
neopulmo  and  to  record  cardiac  output,  ventilation,  and  partial  pressures 
of  inspired  and  expired  gas  and  of  venous  and  arterial  blood.   Comparison 
of  these  recorded  values  with  those  predicted  by  the  model  would  indicate 
the  validity  of  the  model. 

To  determine  possible  receptor  locations,  it  would  be  necessary  to 
first  "calibrate"  the  receptor's  sensitivity  to  carbon  dioxide  partial 
pressure  by  taking  recordings  from  the  vagus  nerve.   Then,  when  the  lung 
is  unidirectionally  ventilated  the  recordings  from  the  vagus  nerve  will 
indicate  the  CO.  partial  pressure  in  the  microenvironment  of  the  receptor 
and  the  model  will  indicate  which  locations  are  compatible  with  that  partial 
pressure. 

Another  possible  experiment,  alluded  to  in  Chapter  VI,  would  be  to 
record  from  a  bird  the  levels  of  cardiac  output  and  ventilation  under 
different  conditions  of  oxygen  uptake.   Since  these  two  variables  can  be 
combined  in  an  infinite  number  of  ways  to  yield  the  same  oxygen  uptake, 
it  would  be  interesting  to  see  which  combination  is  actually  selected  by 
a  particular  bird. 

Future  work  on  the  model  itself  could  include  the  addition  of  a 
compartment,  or  compartments,  to  represent  the  body  tissues.   Within 
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such  a  compartment  the  blood  would  lose  oxygen  and  gain  carbon  dioxide.   In 
this  way  the  partial  pressures  in  the  venous  blood  would  not  have  to  be 
set  arbitrarily,  but  would  be  determined  within  the  model  itself  and  would 
be  dependent  upon  the  constraints  placed  on  the  model. 

The  accuracy  of  Case  5  of  the  model  could  be  improved  by  using  a  third 
order  polynomial  representation  of  the  oxygen  dissociation  curve,  instead 
of  a  second  order  representation  as  has  been  done  in  this  research. 
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APPENDIX  A 

Definitions  of  symbols 

x  location  along  either  capillary  (cm) 
z  location  along  the  parabronchus  (cm) 
1        length  of  capillaries  (cm) 

A        length  of  the  parabronchus  (cm) 

2 
a  cross  sectional  area  of  an  air  capillary  (em  ) 

2 
$        cross  sectional  area  of  the  parabronchus  (cm  ) 

w  capillary  density  along  a  parabronchus  (cm  ) 

P,  partial  pressure  in  the  blood  capillary  (torr) 

P  partial  pressure  in  the  air  capillary  (torr) 

P,  partial  pressure  in  the  n —  cell  of  the  blood  capillary  (torr) 

P  partial  pressure  in  the  n —  cell  of  the  air  capillary  (torr) 

P,  partial  pressure  in  venous  blood  entering  blood  capillary  (torr) 

P.  partial  pressure  in  blood  leaving  blood  capillary  (torr) 

P,  partial  pressure  in  blood  leaving  parabronchus  (torr) 

P  partial  pressure  in  air  capillary  at  open  end  (torr) 

dP  -1 


dx 


derivative  of  P  at  open  end  (torr- cm  ) 


P        partial  pressure  in  the  parabronchus  (torr) 

P  .       partial  pressure  in  air  entering  parabronchus  (torr) 

P        partial  pressure  in  air  leaving  parabronchus  (torr) 

C         concentration.   It  has  the  same  subscripting  as  P  (mmol-lT  ) 

f        rate  of  blood  flow  in  blood  capillary  Q.'im'n  ) 

B        solubility  coefficient  in  blood  capillary  (mmol'L-  •  torr  ) 
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h        diffusing  capacity  between  blood  and  air  capillaries  per  unit 
length  (mmol-min  *torr   -cm  ) 

F  rate  of  gas  flow  in  the  parabronchus  CL-min"  ) 

B  solubility  coefficient  in  parabronchus  (mmol*L   -torr"  ) 

n  diffusion  constant  in  air  capillary  (mmol-min  'cm  'torr  ) 

9  diffusion  constant  in  parabronchus  (mmol-min-  'cm-  -torr"  ) 

Q  total  blood  flow  to  a  parabronchus  (L-min  ) 

Q         flow  of  blood  which  has  not  yet  undergone  exchange  at  same 
location  along  the  parabronchus  (L-min-1) 

Q         flow  of  blood  which  has  undergone  gas  exchange  at  some  location 
along  the  parabronchus  (L»min_1) 
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Publication  Manuscript:   Theory  of  Gas  Exchange  in  the  Avian  Parabronchus. 
William  D.  Crank  and  Richard  R.  Gallagher 


ABSTRACT.   To  understand  the  distribution  of  oxygen  and  carbon  dioxide 
in  the  avian  lung,  a  theoretical  treatment  of  gas  exchange  in  the  para- 
bronchus of  the  avian  lung  is  described.   The  model  used  is  modified 
after  Zeuthen  (1942).   In  addition  to  bulk  flow  through  the.  parabronchial 
lumen,  mass  transport  by  diffusion  through  the  air  spaces  of  both  the 
parabronchial  lumen  and  air  capillaries  is  treated.   The  relationship  of 
the  partial  pressures  of  oxygen  and  carbon  dioxide  within  the  blood 
capillaries,  air  capillaries,  and  parabronchial  lumen  to  parabronchial 
blood  flow  and  ventilation  is  graphically  shown. 

The  results  indicate  that  the  maximum  partial  pressure  variations 
of  oxygen  and  carbon  dioxide  along  the  length  of  an  air  capillary  are 
less  than  one  torr  and  five  torr,  respectively.   Also,  the  resistance 
to  gas  exchange  imposed  by  diffusion  across  the  blood-gas  Interface 
is  many  times  greater  (about  250  times  greater  for  0„;  ranging  from 
8.5  to  33  times  greater  for  CO.)  than  the  resistance  imposed  by  dif- 
fusion through  the  air  capillaries.   The  assumption  of  bulk  flow  only 
through  the  parabronchial  lumen  yields  accurate  results  at  resting  and 
high  ventilation  rates,  but  axial  diffusion  in  the  lumen  must  additionally 
be  considered  at  lower  ventilation  rates.   Before  physically  entering  the 
parabronchus,  the  gas  in  the  pre-parabronchial  region  shows  a  reduction 
in  0.  partial  pressure  due  to  forward  diffusion  toward  the  parabronchus 
and  an  increase  in  CO.  partial  pressure  due  to  backward  diffusion  from 
the  parabronchus.   Also,  the  way  In  which  the  respiratory  exchange  ratio 
varies  with  position  along  the  parabronchus  is  shown. 
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INTRODUCTION 

The  avian  respiratory  system,  in  contrast  to  that  of  mammals  or 
reptiles,  has  rigid  lungs  in  which  gas  is  exchanged  and  air  sacs  which 
act  as  bellows  to  ventilate  those  lungs.   The  parabronchi,  numbering  in 
the  hundreds,  are  the  functional  units  of  each  lung.   Each  parabronchus 
is  basically  a  cylindrical  structure  penetrated  by  a  lumen  along  its 
central  axis.   Ventilated  gas  passes  down  the  lumen,  and  loses  oxygen 
and  picks  up  carbon  dioxide.   Oxygen  diffuses  into  the  air  capillaries, 
which  radiate  outward  from  the  parabronchial  lumen,  and  carbon  dioxide 
diffuses  in  the  opposite  direction.   Within  an  air  capillary  there  is  no 
bulk  flow  of  gas  (Zeuthen,  1942).   Gas  is  exchanged  across  the  surfaces 
separating  air  capillaries  from  the  adjacent  blood  capillaries.   The 
direction  of  flow  through  the  blood  capillaries  is  from  the  parabronchial 
periphery  inward  toward  the  parabronchial  lumen  (Abdalla  and  King,  1974; 
West,  Bamford,  and  Jones,  1977).   Venous  blood  enters  the  capillaries 
at  the  periphery  and  the  blood  leaving  the  capillaries  mixes  to  form 
arterial  blood. 

Scheid  and  Piiper  (1970)  developed  a  cross-current  model  to  describe 
the  function  of  the  avian  parabronchus.   In  their  model,  gas  is  exchanged 
between  the  gas  moving  by  bulk  flow  through  a  tube  and  a  thin  sheet  of 
blood  flowing  across  the  surface  of  that  tube.   The  direction  of  blood 
flow  is  perpendicular  to  the  direction  of  gas  flow,  which  is  the  reason 
they  describe  the  model  as  "cross-current." 

In  the  cross-current  model  it  is  assumed  that  the  partial  pressure 
of  O2  or  CO2  is  constant  throughout  a  cross-section  cutting  the  tube 
perpendicular  to  its  axis.  Applied  to  the  avian  parabronchus,  this 


assumption  forces  the  partial  pressure  along  the  radially  projecting 
air  capillaries  to  be  constant,  which  implies  zero  diffusion  resistance. 
It  is  also  assumed  that  gas  moves  through  the  parabronchial  lumen  only  by 
bulk  flow  without  axial  diffusion,  which  implies  infinite  diffusion 
resistance.   Diffusion  through  the  air  spaces  must  occur,  however,  with 
equal  resistance  in  both  the  radial  and  axial  directions  if  a  square 
velocity  profile  across  the  lumen  is  assumed. 

In  the  model  described  here,  the  blood  and  air  capillaries  are 
retained  as  separate  structures.   Also,  the  diffusion  constant  of  0.  or 
or  CO.  in  air  characterizes  the  diffusion  radially  along  the  air  capil- 
laries and  axially  along  the  lumen.   Constant  partial  pressure  is  still 
assumed  for  a  cross  section  through  the  lumen. 

THEORY 

The  physical  structure  of  the  model  is  shown  in  Figure  1A.  This 
representation  is  similar  to  that  of  Zeuthen  (1942).   Explicitly,  the 
more  important  assumptions  of  the  model  are:   1)  Steady  state  operation, 
2)  Homogeneous  perfusion  of  the  parabronchus ,  3)  Homogenous  exchange 
surface  between  gas  and  blood,  4)  Linear  0.  and  CO.  dissociation  curves, 
5)  No  interaction  between  0,  and  CO.,  6)  Square  velocity  profile  across 
the  parabronchial  lumen. 

All  the  blood  which  exchanges  gas  with  a  given  air  capillary  is 
taken  to  flow  through  a  single,  adjacent  blood  capillary.   It  is  assumed 
that  the  air  capillaries  are  the  same  length  as  the  blood  capillaries. 
Both  air  and  blood  capillaries  are  represented  as  a  series  of 
homogeneous  compartments  of  length  Ax  with  mass  transport  across  com- 
partment surfaces.   The  expressions  describing  this  transport  between 


VENOUS  BLOOD 

AIR  CAPILLARY 1 

BLOOD  CAPILLARY 1 

PARABRONCHIAL 
LUMEN  VI    I 


compartment 


Figure  1.   A:   The  parabronchial  model,  modified  after  Zeuthen.   B: 

Detail  of  air  capillary  and  blood  capillary  compartments.   Expressions 
describing  transport  of  the  considered  gas  between  compartments  are 
shown.   C:   Detail  of  parabronchial  luminal  compartments.   Expressions 
are  shown  which  describe  transport  between  luminal  compartments  by 
both  bulk  flow  and  diffusion  and  transport  between  a  luminal  compartment 
and  the  air  capillaries  which  open  into  it. 
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compartment  n  and  the  adjacent  compartments  are  shown  in  Figure  IB.   The 
symbols  used  are  defined  in  Appendix  C.   In  a  single  air  or  blood 
capillary  compartment,  equating  the  influx  of  the  considered  gas  to  the 
efflux  yields  a  difference  equation.   Performing  the  limit  Ax  +  0  trans- 
forms this  into  a  differential  equation.   In  Appendix  A  the  derivation 
of  differential  equations  for  air  capillary  and  blood  capillary  partial 
pressures  and  the  boundary  conditions  on  the  solutions  of  those  equations 
are  given  in  detail. 

Illustrated  in  Figure  1A  are  three  regions,  denoted  I,  II,  and  III, 
through  which  the  parabronchial  gas  consecutively  flows.   The  parabronchus 
is  coincident  with  region  II.   The  parabronchial  entrance  and  exit  lie 
at  z=0  and  z=L.   Although  this  geometry  only  approximately  depicts  the 
opening  of  a  parabronchus  into  a  secondary  bronchus,  the  results  should 
give  some  insight  into  the  gas  transport  processes  occurring  there. 

By  a  compartmentalization  similar  to  that  of  the  air  and  blood 
capillaries,  the  parabronchial  lumen  is  considered  to  be  composed  of  a 
series  of  homogeneous  compartments  of  length  Az.   Transport  between 
compartment  n,  adjacent  compartments,  and  the  air  capillaries  is  illus- 
trated in  Figure  1C.   The  transport  due  to  bulk  flow  is  shown  separately 
from  that  due  to  diffusion.   The  transport  due  to  all  the  air  capillaries 
opening  into  the  compartment  accounts  for  the  gas  exchanged  between  the 
blood  phase  and  gas  phase.   In  a  single  parabronchial  compartment,  equating 
the  influx  of  the  considered  gas  to  the  efflux  and  performing  the  limit 
Az  -*■  0  yields  a  differential  equation  for  partial  pressure  along  the 
parabronchus.   The  differential  equations  for  partial  pressure  in  regions 
I,  II,  III  and  the  boundary  conditions  on  their  solutions  are  given  in 
Appendix  A. 
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In  Appendix  B  the  differential  equations  are  solved  for  the  partial 
pressure  in  the  air  capillaries,  blood  capillaries  and  regions  I,  II,  III 
along  the  parabronchial  lumen. 
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air  capillary  have  been  calculated  from  equations  (29)  and  (30)  of 
Appendix  B.  Figure  2A  illustrates  that  within  a  blood  capillary  the 
0,  partial  pressure  increases  from  the  venous  end  at  the  parabronchial 
periphery  (X=0)  to  the  arterial  end  at  the  parabronchial  lumen  (X=0.024). 
Figure  2B  illustrates  that  within  an  air  capillary  the  0.  partial  pressure 
increases  from  the  closed  end  at  the  periphery  (X=0)  to  the  open  end  at 


pressure  from  the  venous  end  to  the  arterial  end  of  the  blood  capillary. 
Figure  2D  illustrates  the  decrease  of  CO,  partial  pressure  from  the 
peripheral  end  to  the  luminal  end  of  the  air  capillary.   This  gradient 
forces  CO  to  diffuse  from  the  air  capillary  into  the  lumen. 

In  Figure  2  the  air  capillary  has  been  assumed  to  open  into  para- 
bronchial gas  with  P  .=160  and  Prn2=0  torr.   Two  cases  have  been  con- 
sidered for  the  parabronchial  perfusion:   Q  =0.29  and  Q  =2.9  ml/min.   Two 


P   =20  and  P   =40  torr.   Likewise,  two  cases  for  the  CO.  partial  pressure 

V02         V02  l 

of  the  venous  blood  have  been  taken  as:   P    =45  and  P    =55  torr. 

VC02         VC02 
For  both  0,  and  CO.,  in  Figures  2A  and  2C,  increasing  the  capillary 

blood  flow  decreases  the  partial  pressure  variation  along  the  length  of 

the  blood  capillary.   But  for  both  gases,  in  Figures  2B  and  2D,  increasing 

blood  flow  increases  the  partial  pressure  variation  along  the  length  of 

the  air  capillary.   An  increased  partial  pressure  gradient  in  the  air 

capillary  means  that  mass  transport  by  diffusion  is  increased.   So,  when 


Figure  2.   Partial  pressures  along  a  blood  capillary  and  an  air  capillary 

when  Q  =0.29  ml/min  (dotted  lines)  and  Q  =2.9  ml/min  (solid  lines). 

P  P 

The  cases  considered  for  venous  partial  pressure  are  P-Q  =20,  P-Q  =40, 

p-   =45,  and  P-   =55  torr.   The  luminal  gas  at  the  open  end  of  the 

air  capillary  is  described  by  PQ  =160  and  PCQ  =0  torr.   A:   PQ  along 

the  blood  capillary.   B:   P   along  the  air  capillary.   C:   PCQ  along 

the  blood  capillary.   D:   P    along  the  air  capillary. 
LU2 


the  capillary  blood  flow  increases,  the  exchange  of  0.  and  CO.  between 

the  blood  phase  and  the  parabronchial  gas  increases,  even  though  the 

blood  equilibration  decreases. 

When  P    decreases  without  a  change  in  capillary  blood  flow,  a 
V02 
greater  amoung  of  0.  is  moved  from  the  parabronchial  gas  into  the  blood. 

In  Figures  2A  and  2B  this  is  shown  by  an  increase  in  the  partial  pressure 

variation  along  both  the  air  capillary  and  blood  capillary.   Likewise, 

Figures  2C  and  2D  show  that  when  P     increases  without  a  change  in 

VC02 
capillary  blood  flow,  the  partial  pressure  variation  increases  along  both 

the  air  capillary  and  blood  capillary. 

Comparison  of  Figures  2A  and  2C  shows  that  the  capillary  blood  CO, 


In  alveolar  lungs  the  rate  at  which  pulmonary  capillary  blood  equilibrates 
with  alveolar  gas  depends  upon  the  ratio  of  the  diffusing  capacity  to 
the  capacitance  coefficient  of  blood  (Forster,  1964;  Wagner,  1977).  This 
same  ratio  (see  the  definition  of  p  in  Appendices  B  and  D)  enters  into 
the  equations  of  Appendix  B  which  describe  the  solution  for  partial 


coefficient  of  CO.  in  blood  is  taken  larger  than  that  of  0,  by  a  factor  of 
3.36,  the  ratio  of  these  quantities  is  larger  for  CO.  than  for  0.  by  a 
factor  of  5.95.  Thus,  the  rate  of  exponential  equilibration  along  the 
blood  capillary  is  greater  for  CO.  than  for  0. . 

The  variation  in  0.  partial  pressure  along  the  length  of  the  air 


partial  pressure  is  less  than  five  torr  (Figure  2D).   For  both  0.  and 
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CO.,  the  high  Q  values  and  extreme  parabronchial  partial  pressures  have 

been  chosen  to  give  partial  pressure  variations  along  the  air  capillary 

which  are  near  the  maximum  of  the  normal  physiological  range.   The 

capillaries  which  experience  conditions  most  nearly  approaching  these 

lie  at  the  parabronchial  entrance.   The  following  capillaries  along  the 

parabronchus  experience  luminal  gas  which  has  already  undergone  some 

exchange.   So  the  driving  force  for  further  gas  exchange  is  reduced  and 

the  partial  pressure  variation  along  these  air  capillaries  is  less  than 

that  shown  in  2B  and  2D. 

Figures  3A  and  3B  show  that  along  the  parabronchus  from  the  entrance 

to  the  exit,  the  0     partial  pressure  decreases  and  the  CO.  partial  pressure 

increases.   There  are  two  families  of  curves  in  both  3A  and  3B.   The  family 

of  dotted  curves  was  computed  from  equations  (47),  (48),  and  (49)  of 

Appendix  B  under  the  assumption  of  equal  diffusion  resistance  in  the  air 

capillaries  and  the  parabronchial  lumen.   The  family  of  solid  curves  was 

computed  under  the  assumption  of  bulk  flow  only  in  the  parabronchial 

lumen  (infinite  diffusion  resistance)  and  constant  partial  pressure  along 

an  air  capillary  (zero  diffusion  resistance) . 

In  Figure  3  the  parabronchial  perfusion  is  Q  =0.29  ml/min  and  the 

P 

parabronchial  ventilation  varies  from  V  =0.022  to  V  =0.90  ml/min.   The 

P  P 

state  of  venous  blood  is  given  by  P   =40  and  P    =45  torr.   The  state 

V02         VC02 
of  inspired  gas  is  given  by  P   =160  and  P    =0  torr.   The  diffusion 

constants  are  K  .,=6.28  x  10~  and  K  .=4.90  x  10~  mmol-min"  -torr"  -cm"  . 

The  dotted  curves,  in  which  diffusion  is  accounted  for,  show  a 

change  in  the  partial  pressure  of  the  ventilating  gas  before  it  enters 

the  parabronchus  (region  II) .   This  is  in  contrast  to  the  solid  curves 
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Figure  3.   Partial  pressures  along  the  parabronchus  computed  under  the 

assumption  of  (1):   equal  diffusion  resistance  in  the  air  capillaries 
and  parabronchial  lumen  (dotted  curves)  and,  (2):   constant  partial 
pressure  along  the  air  capillaries  and  bulk  flow  only  in  the  para- 
bronchial lumen  (solid  curves).   The  parabronchial  ventilation  varies 


i02 
=160  and  P,rf)  =0  torr;  the  state  of  venous  blood  is  V-     =40  and 

P-r  *=45  torr.   A:   P   along  the  parabronchus.   B:   P    along  the 

parabronchus . 


B12 


which  neglect  diffusion.   Figure  3A  shows  that  the  partial  pressure  in 
the  pre-parabronchial  region  decreases  because  0.  diffuses  out  of  this 
region  into  the  parabronchus  which  has  a  lower  00  partial  pressure. 
Figure  3B  shows  that  the  partial  pressure  in  the  pre-parabronchial  region 
increases  because  CO.  diffuses  backward  into  this  region  from  the  para- 
bronchus which  has  a  higher  CO.  partial  pressure. 

Figures  3A  and  3B  show  that  increasing  the  parabronchial  ventilation 
causes  the  gas  to  equilibrate  less  fully  with  venous  blood.   At  high 
ventilation  rates  the  curves  computed  with  parabronchial  diffusion  are 
closely  approximated  by  the  curves  computed  with  bulk  flow  only.   At  lower 
ventilation  rates  this  approximation  no  longer  holds  as  diffusion  becomes 
relatively  more  important  in  comparison  with  bulk  flow. 

The  partial  pressures  computed  with  parabronchial  diffusion  show  less 
change  between  the  inspired  and  end-parabronchial  values  than  do  the 
partial  pressures  computed  with  bulk  flow  only.   Depending  on  the  ventila- 
tion rate,  that  change  is  less  by  factors  as  small  as  0.94  and  0.98  for 
oxygen  and  carbon  dioxide,  respectively. 

Comparison  of  3A  and  3B  shows  that  at  a  fixed  ventilation  the  rate 
of  equilibration  along  the  parabronchus  is  greater  for  CO.  than  for  0_. 
This  is  a  secondary  consequence  of  the  characteristics  of  gas  exchange 
along  a  blood-air  capillary  pair.   As  shown  in  Figures  2B  and  2D,  at 
the  parabronchial  entrance  the  partial  pressure  gradient  along  the  air 


Figure  4).   This  further  implies  that  in  the  parabronchial  lumen,  CO. 
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Figure  4  gives  a  more  complete  picture  of  the  whole  parabronchus 
as  a  gas  exchange  organ  by  combining  the  results  describing  partial 
pressure  variation  along  the  blood  capillaries  with  those  describing  the 
variation  along  the  parabronchial  lumen.   Results  for  0-  are  shown  in  4A; 
for  C0~  in  4B.  The  series  of  curves  on  the  left  describes  partial  pres- 
sure variation  in  blood  capillaries  located  at  different  positions  along 
the  parabronchus.   In  each  capillary  the  partial  pressure  changes  from 
the  venous  value  to  the  end-capillary  value.   The  rate  and  extent  of  that 
change  depends  on  the  partial  pressure  of  the  gas  in  the  adjacent  air 
capillary.   The  partial  pressure  in  the  air  capillary  depends  upon  the 
partial  pressure  in  the  parabronchial  gas  into  which  the  air  capillary 
opens,  which  in  turn  depends  upon  its  location  along  the  parabronchus. 
The  variation  in  partial  pressure  of  the  gas  along  the  parabronchus  is 
shown  by  the  curve  on  the  right.   The  physical  location  of  a  blood  capil- 
lary with  respect  to  position  along  the  parabronchus  is  indicated  by  the 
dotted  lines.   The  vertical  segment  of  the  dotted  line  indicates  the  extent 
to  which  the  capillary  blood  falls  short  of  full  equilibration  with  the 
gas  in  the  adjacent  air  capillary.   Gas  leaving  the  parabronchus  is  at 
the  partial  pressure  indicated  by  the  right  end  point  of  the  parabronchial 
gas  curve.   The  partial  pressure  of  arterial  blood,  formed  from  the  blood 
emptying  out  of  all  the  blood  capillaries  along  the  parabronchus,  is  an 
average  of  the  values  of  partial  pressure  of  the  right  end  points  of  the 
blood  capillary  curves. 

Figure  5  has  been  constructed  in  order  to  compare  the  rate  of  0 
uptake  (MQ2)  to  the  rate  of  C02  elimination  (Mc0,)  by  a  single  blood 
capillary.   The  computed  quantities  MQ2,  M   ,  and  ^CQ2^Q2    (che 
respiratory  exchange  ratio)  have  been  plotted  against  distance  along  the 
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Figure  4.   A:   02  partial  pressure  along  several  parabronchial  blood 
capillaries  (left)  and  along  the  parabronchial  lumen  (right).   The 
capillary  location  along  the  parabronchus  is  indicated  by  horizontal 
dotted  lines.   The  extent  to  which  blood  falls  short  of  full  equili- 
bration with  air  is  indicated  by  vertical  dotted  lines.   B:   CO 

2 

partial  pressure  along  several  parabronchial  blood  capillaries  and 
along  the  parabronchial  lumen. 
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Figure  5.   Gas  exchange  rate  of  CO  and  0  by  a  single  blood-air  capillary 
pair  plotted  against  distance  along  the  parabronchus.   Also  shown  is 
the  ratio  of  these  exchange  rates. 
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parabronchus.   In  order  to  quantify  the  single  capillary  gas  exchange 

rate,  the  number  of  blood-air  capillary  pairs  per  parabronchus  has  been 

taken  to  be  3.1  x  10  . 

Since  the  capillaries  near  the  parabronchial  entrance  experience 

the  freshest  gas  and  the  following  capillaries  experience  luminal  gas 

which  has  already  under 

along  the  parabronchus. 

*  b 
C02 

causes  the  respiratory  exchange  ratio,  M  _/M  -,  also  to  decrease  along 
the  parabronchus.   It  changes  from  1.3  for  the  capillary  pairs  at  the 
entrance  to  0.21  for  those  at  the  exit.   The  over-all  respiratory  exchange 
ratio  for  this  entire  parabronchus  is  0. 61.   This  value  is  less  than  the 
expected  value  of  about  0.8,  probably  because  the  model's  use  of  a  linear 
0?  dissociation  curve  exaggerates  the  0-  exchange. 
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DISCUSSION 

The  assumption  of  a  linear  dissociation  curve  is  an  approximation 
for  the  respiratory  gases,  oxygen  and  carbon  dioxide.   A  consequence  of 
the  linear  approximation  for  the  sigmoid  09  dissociation  curve  over  the 
normal  physiological  range  is  an  unrealistically  high  0~   concentration 


Therefore,  the  linear  dissociation  curve  used  here  results  in  an  exag- 
geration of  the  computed  0_  partial  pres 
capillaries  and  along  the  parabronchus. 

In  contrast  to  0-,  the  effective  CC 
approximated  by  a  straight  line  (Schiepers,  Kawashiro,  and  Scheid,  1975). 

A  square  velocity  profile  across  the  parabronchial  lumen  has  been 
assumed  in  the  calculations  of  partial  pressure  along  the  parabronchus. 
If  the  actual  velocity  profile  deviates  from  this  ideal,  there  will  be 
additional  mixing  of  gases  in  the  lumen  which  is  not  accounted  for. 
This  implies  an  under-estimation  of  the  effective  diffusion  in  the 
curves  of  Figure  4. 

Blood  flows  in  the  capillaries  from  the  periphery  of  the  parabronchus 
toward  the  lumen.   In  the  air  capillaries  lying  parallel  to  the  blood 
capillaries,  oxygen  diffuses  from  the  lumen  toward  the  periphery  and 
carbon  dioxide  diffuses  from  the  periphery  toward  the  lumen.   This  means 
that  in  the  blood-air  capillary  network  there  are  counter-current  oxygen 
flows  occurring  simultaneously  with  co-current  carbon  dioxide  flows. 
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In  the  parabronchus,  the  radial  flow  of  capillary  blood  is  perpen- 
dicular to  the  axial  flow  of  luminal  gas.   Thus,  the  over-all  parabronchial 
process  can  be  described  as  cross-current  gas  exchange.   The  0_  and  CO, 
diffusive  flows  in  the  air  capillaries,  which  are  counter-current  and 
co-current  with  respect  to  the  0.  and  CO,,  bulk  flows  in  the  blood 
capillaries,  move  these  gases  between  the  cross-current  flows  of 
luminal  gas  and  capillary  blood. 

Between  the  capillary  blood  and  the  parabronchial  luminal  gas  there 
is  a  resistance  to  diffusion  across  the  interface  separating  the  blood 
phase  from  the  gas  phase  and  there  is  a  resistance  to  diffusion  through 
the  gas  phase  of  the  air  capillaries.   We  want  to  obtain  an  estimate  of 
the  relative  importance  of  interface  resistance  to  air  capillary  resis- 
tance within  a  single  blood-air  capillary  pair.   This  estimate  is  obtained 

by  evaluating  the  expression:   (P,    -P    )/(P    -P__).   P,    is  an 

b,av  g,av    g,av  II     b,av 

average  partial  pressure  in  the  blood  capillary;  P     is  an  average  partial 

pressure  in  the  air  capillary;  P   is  the  partial  pressure  at  the  luminal  end 

of  the  air  capillary.   Partial  pressures  are  taken  as  average  at  that 

point  along  the  blood-air  capillary  pair  where  the  gas  exchange  is  half 

completed.  The  rationale  behind  this  expression  is  that  at  a  given  rate 

of  diffusive  flow,  the  diffusion  resistance  between  two  points  is  propor- 

tonal  to  their  partial  pressure  difference. 

Partial  pressure  values  from  the  curves  in  Figure  2  have  been  used 

to  evaluate  the  above  expression.   From  2A  and  2B  it  is  found  that  the 

resistance  to  0,  uptake  imposed  by  diffusion  from  the  air  capillary  into 

the  blood  capillary  is  about  250  times  the  resistance  imposed  by  diffusion 

through  the  air  capillary  for  the  four  conditions:   1)  P_  =20  torr, 

V02 
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Q  =0.29  ml/mlnj  2)  P   =40  torr,  Q  =0.29  ml/min;  3)  P   =20  torr,  Q  =2.9 

P  V02   _        p  V02  p 

ml/min;  4)  P_  =40  torr,  Q  =2.9  ml/min.   From  2C  and  2D  it  is  found  that 


V02 


the  resistance  to  CO.  excretion  imposed  by  diffusion  from  the  blood 

capillary  into  the  air  capillary  is  about  33  times  the  resistance  imposed 

by  diffusion  through  the  air  capillary  for  the  two  conditions:   1)  P    =45 

VC02 
torr,  Q  =0.29  ml/min;  2)  P    =55  torr,  Q  =0.29  ml/min.   It  is  found  to 

P  VC02  p 

be  about  8.5  for  the  other  two  conditions:   3)  P    =45  torr,  Q  =2.9  ml/min; 

VC02  p 

4)  P    =55  torr,  Q  =2.9  ml/min. 
VC02  P 

In  contrast  to  the  conclusions  of  Piiper  and  Scheid,  1977,  these 

results  show  a  difference  between  0_  and  CO.  in  the  comparison  of  the 

importance  of  interface  diffusion  resistance  to  air  capillary  diffusion 

resistance  to  limit  gas  exchange.   The  principal  reason  for  the  difference 

between  the  two  respiratory  gases  in  this  comparison  is  that  the  diffusing 

capacity  of  CO.  across  the  blood-gas  interface  is  20  times  greater  than 

that  of  0.. 

To  find  the  relative  importance  of  interface  resistance  to  air 
capillary  resistance  for  the  whole  parabronchus,  it  would  be  necessary 
to  average  the  effects  of  all  the  blood-air  capillary  pairs  along  the 
entire  length  of  the  parabronchus. 

The  results  shown  in  Figure  3  predict  that  the  gas  in  the  pre- 
parabronchial  region  will  begin  to  show  a  reduction  in  0.  partial  pres- 
sure and  an  increase  in  CO^  partial  pressure  before  entering  the  para- 
bronchus.  This  is  due  to  forward  diffusion  of  0.  toward  the  parabronchus 
and  backward  diffusion  of  CO,  from  the  parabronchus.   Measurements  of 
C02  partial  pressure  from  the  dorsobronchi  of  ducks,  in  which  the 
ventilation  is  unidirectional  during  the  respiratory  cycle  but  variable 
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in  flow  rate  (see  for  instance  Scheid  and  Piiper,  1971),  have  been  made 
by  Jammes  and  Bouverot  (1975) .   These  measurements  show  phasic  changes 
during  the  respiratory  cycle,  but  the  authors  attribute  these  entirely 
to  bulk  movement  of  gas.   Additional  measurements  will  be  necessary  to 
determine  the  importance  of  backward  diffusion  on  dorsobronchial  CO. 
partial  pressure  during  respiration  in  a  live  bird. 

An  assumption  made  by  Scheid  and  Piiper  (1970)  in  the  analysis  of 
their  cross-current  model  was  that  the  partial  pressures  of  0_  and  CO. 
are  constant  throughout  the  gas  phase  of  a  cross  section  cutting  the 
model  parabronchus  perpendicular  to  its  axis.   A  second  assumption  made 
was  the  absence  of  any  transport  in  the  parabronchial  lumen  by  diffusion 
in  the  axial  direction.   In  the  parabronchial  model  shown  in  Figure  1  the 
conditions  of  the  first  assumption  are  reached  in  the  limit  of  an  infinitely 
large  diffusion  constant  for  0,  and  CO,  within  the  air  capillaries.   The 
conditions  of  the  second  assumption  are  reached  in  the  limit  of  a  zero 
valued  diffusion  constant  for  0,  and  CO,  within  the  parabronchial  lumen. 
Under  these  two  limiting  processes  the  solution  for  P   (z)  (equation 
(48)  in  Appendix  B) ,  when  evaluated  at  the  parabronchial  exit,  reduces 
to  equation  (7)  of  Scheid  and  Piiper. 

The  results  shown  in  Figure  2  support  the  reasonableness  of  the 
first  simplifying  assumption  except  for  CO,  in  the  neighborhood  of  the 
parabronchial  entrance  under  conditions  of  high  blood  flow  rate.   The 
results  shown  in  Figure  3  support  the  second  simplifying  assumption 
except  under  conditions  of  low  ventilation  rate. 
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APPENDIX  A.   The  Differential  Equations 

Under  the  steady-state  conditions  in  air  capillary  compartment  n, 
the  influx  of  the  considered  gas  equals  its  efflux.   This  fact,  together 
with  the  expressions  of  Figure  IB,  leads  to  the  equation 

(1)  fK[?^n+^  "  PS>"  -  PS,"  "  V"-!]  .  h(P    .  p   )4x  . 

Ax  Ax  g,n    b,n 

Taking  the  limit  Ax-K)  in  (1)  yields  the  differential  equation 

d2P 

(2)  fK  |  =  h(W  • 

dx 

Similar  considerations  applied  to  blood  capillary  compartment  n  lead  to 
the  differential  equation 

dP. 

(3)  Sq  2J*  =  h(Pg-Pb)  . 

Equating  the  left  members  of  (2)  and  (3)  yields 

d2P       dP 

(4)  fK  —a  -  Bq  ^  • 

dx 

The  boundary  conditions  on  solutions  of  (2),  (3),  and  (4)  are  imposed 
at  the  ends  of  the  capillary  pair.   Both  the  closed  end  of  the  air 
capillary  and  the  point  at  which  venous  blood  enters  the  blood  capillary 
lie  on  the  parabronchial  periphery  at  x=0.   The  open  end  of  the  air 
capillary  and  the  point  at  which  blood  leaves  the  blood  capillary  lie 
on  the  parabronchial  luminal  surface  at  x=&.   The  boundary  conditions 


(5) 


at  x=0: 

Pb(0) 

=  pv 

at  x=0: 

dP    , 

& 

dx 

=       ( 

(6)  at  x=0:         -^     ■  0  , 

~'x=0 
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(7)  at  x=£: 

Boundary  condition  (6)  results  from  the  lack  of  diffusion  across  the  end 

of  the  air  capillary.   Integration  of  (4)  over  the  interval  [0,x]  and 

application  of  the  boundary  conditions  (5)  and  (6)  yields 

dP 
(8) 

Setting  x=£  in  (8),  one  obtains 

dP  I 

(9)  fK  — &     =  gq(P  OD-P  )  . 

'x=l 

Either  the  right  or  left  member  of  (9)  describes  the  transport  between  a 
single  air  capillary  and  the  parabronchus.   The  transport  due  to  all  the 
air  capillaries  opening  into  a  single  parabronchial  compartment  is 
obtained  from  the  right  hand  side  of  (9)  by  replacing  q  with  (Q  Az/L) . 
Assuming  steady-state  conditions  for  parabronchial  compartment  n 
(Figure  1C),  equating  the  influx  to  the  efflux,  and  taking  the  limit 
Az-+0,  leads  to  the  differential  equation 

(10)  in  region  II:     FK  ^j  -  aV^  %  =  -3-   (P  00  -  P..)  . 

.  L  p  dz    L    b       V 

dz 

In  regions  I  and  III  there  is  no  gas  exchange,  so  the  differential 

equation  to  be  satisfied  is 

2 

(11)  in  regions   I  and   III:      FK  =-|  -  ctV     ^-  =   0    . 

A,2  P   dz 


(12) 
(13) 


(14)  at   z=0    :  _£  -  -JU 

dz  dz 
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(15) 
(16) 
(17) 


II 

dP 


II 
dz 


III    ' 

dPIII 
dz 


iV-r 


=   0 
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APPENDIX  B.   Solutions  of  Differential  Equations 

Next,  the  solutions  for  partial  pressure  along  the  blood  and  air 

capillaries  are  found.   Taking  the  Laplace  transforms  of  equations  (3) 

and  (8) ,  one  obtains  two  equations  in  the  two  unknowns  P  and  P,  .  which 

g      b 


(18) 

8        s     +  sp 

—  y°>  -     2 

-  a                   s(s     +  sp   - 

-  a) 

(19) 

P       i                 p 

2 

r    mi    i           s     -  a 

s       +    Sp 

r  (U)   i        2 
-  a                    s(s     +  sp   - 

-  a) 

where 

(20) 

h 

a                h 
and       a  =  —  . 

Taking  inverse  transforms  of  (18)  and  (19)  yields  solutions  for  P  and 

V 

(21)  P  (x)  =  H(x)  P  (0)  +  B(x)Pv  . 

(22)  Pb(x)  =  G(x)  P  (0)  +  D(x)Py  . 

The  functions  H(x),  B(x),  G(x),  and  D(x)  are  given  as  follows. 

S-.  +  p     3  X    S,  +  p     S.X 

(23)  H(x)  ■  -i e  L     +  -= e  l        , 

Sl   S2        s2  ~  Sl 

s,x  s^x 

(24)  B(x)  =  1  - 


(25)      G(x) 


(26)      D(x)  -  1  + 


s1(s1-s2) 

S.  X 
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where 

(27)      sx  -  =|  -  |  (p2  +  4a)1/2   ,   and   s2  =  =^  +  \   (p2  +  4a)1/2  . 

A  difficulty  yet  remaining  with  the  solutions  (21)  and  (22)  is  that  they 


set  equal  to  I   in  equation  (21) .   Then  solving  for  P  (0)  yields 

(28)  y°>  -  ik  pu  -  W>  pv  • 

Substituting  (28)  into  (21)  and  (22)  gives  the  solutions  for  P  and  P. 

g      b 

in  the  form 

(29)  P    (x)    =  5$4  P        +    fB(x)    -  H(x)BWlP 
KZV}  VgW         H(z)    Pjj   +    IBW  H(i)     JPV    , 

(30)  P    (x)   =  §££  p        +    rDfx)    -  GOOBO). 

Equations  (29)  and  (30)  are  the  desired  solutions  for  partial  pressure 
along  the  air  and  blood  capillaries. 

Next,  the  solutions  for  partial  pressure  along  the  parabronchus  are 
found.   From  (30)  one  obtains 

(3D  P  U)  -  P    =  ^-p       i    fp(i)      G<*>B<*>    .  nr 

ujj  fbu;      fv     nm  fxi  +  idu;  H(£)    -  ijpv 

Substitution  of  (31)  into  (10)  yields  the  differential  equation  which  must 
be  satisfied  in  region  II. 

(32)  4-f-*P-«v. 
dz 

The  following  substitutions  have  been  made. 

aV 

(33)  t  =  — 2-   , 
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(34)  *=-±Wl 

K3   '  V       FKL  H(4) 


(35)  6   =  — E-  [D(l)    -  SffliW    .   ,, 
The  solution  of    (32)   is 

tiz  c?z        X 

(36)  Pn(z)   =  bxe         +b2e2     -  ±  ^  , 


,„>        ,,      _   t        1.2..  ,.1/2  .      _  x    .    1    ,   2    ,     .  ,.1/2 

(37)  t.^  =  y  -  j  (t     +  4i|i)  ,      and      t2  =  i  +  J  (T     +  **) 

The  equation  which  must   be  satisfied   in  regions   I  and   III   is 

(38)  &  -  r  f  =  0    , 
dz 

where  (33)  has  been  combined  with  (11)  to  arrive  at  this  form.   The 


(39) 

(40) 

The  constants  b  -b,  must  be  chosen  so  that  the  boundary  conditions  are 
satisfied.   Substituting  the  solutions  (36),  (39),  and  (40)  into  boundary 
conditions  (12)-(17)  yields  the  following  equations. 

(41)  b4  =  P.  , 

(42)  b3  +  b,  -  bl  +  b,  -  i  Pv  , 
(43) 


(44)  b^  ♦  b2et2L  -  |  Py  -  b5eTL  +  b,  , 

(45)  tjbje  L     +  t2b2e  *  -  Tb5eTL  , 

(46)  b5  =  0  . 

Solving  these  equations  for  the  b's  and  substituting  the  resulting  values 
into  (36),  (39),  and  (40)  produces  the  desired  solutions  for  P  ,  PTT,  and 


t  L    t  L       , 
t  t  eTZ(e    -  e   )  (P  +  I  P  ) 

(47)  Pi(z)=pi  +  -AA__ .  J-   » 

txe        Ctj-T)  -  t2e   (tj-f) 

t.z+t-L      t.L+t,z 
t(t  e        -  e.«  X     )(P.  +  £pv) 

(48)  P   (z)  =  —  P  +  - *     V 

W8'     tIItZ;    «  ^V  +         t  L  t  L 

'l6    (t2"T)  "  t2e        <-ti~T'> 

L(t  +t  ) 
.      T(t,-t.)e       (P.  +  f  P  ) 

(49)  P   (z)  -  —  P  + — i ^  v 

t1e    (t2-x)  -  t2e    (tx-T) 
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APPENDIX  C  -  Symbol  Definitions 

x        location  along  capillary  (cm) 

I  length  of  capillaries  (cm) 

2 
f        cross-sectional  area  of  air  capillary  (cm  ) 

z         location  along  parabronchus  (cm) 

L        length  of  parabronchus  (cm) 

2 
F         cross-sectional  area  of  parabronchial  lumen  (cm  ) 

h         diffusing  capacity  between  blood  and  air  capillaries  per  unit 
length  (mmol-min  'torr  'cm  ) 

K  Krogh's  diffusion  constant  in  gas  phase  (mmol-min  'torr  'cm  ) 

S  capacitance  coefficient  in  blood  (mmol-ml-  -torr-  ) 

a  capacitance  coefficient  in  gas  phase  (mmol-ml   -torr   ) 

q  blood  supply  to  single  blood  capillary  (ml-min~  ) 

Q  blood  supply  to  single  parabronchus  (ml-min-  ) 

V 

P 

P  partial  pressure  in  cell  n  of  air  capillary  (torr) 

P.  partial  pressure  in  cell  n  of  blood  capillary  (torr) 

P  partial  pressure  as  a  function  of  x  in  air  capillary  (torr) 

Pu  partial  pressure  as  a  function  of  x  in  blood  capillary  (torr) 

PI'PII'I'lII    partial  pressure  as  a  function  of  z  in  region  I,  II,  III 
(torr) 

P^  partial  pressure  of  inspired  gas  (torr) 

P„  partial  pressure  of  venous  blood  (torr) 

s  Laplace  transform  variable 

D  diffusing  capacity  of  single  parabronchus  (mmol-min-  •torr"  ) 

A        cross  sectional  area  of  all  air  capillaries  in  a  single  parabronchus 
P  9 
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APPENDIX  D.  Parameter  Values 

To  obtain  the  parameter  values  used  in  evaluating  the  model  results, 
the  various  quantities  from  different  sources  have  been  adjusted  to  fit 
a  1.6  kg  bird. 

King  and  Payne  (1964)  report  the  minute  ventilation  of  resting 
chickens  (3.4  kg,  female)  to  be  766  ml/min.   Vogel  and  Sturkie  (1963) 
report  the  cardiac  output  of  resting  chickens  (1.95  kg,  female)  to  range 
from  234  ml/min  (summer)  to  345  ml/min  (winter).   From  these  data  repre- 
sentative rates  of  ventilation  and  perfusion  for  a  single  parabronchus 

are  found  to  be  V  =0.45  ml/min  and  Q  =0.29  ml/min  when  the  total  number 
P  P 

of  parabronchi  in  both  lungs  is  taken  to  be  800  (King  and  Cowie,  1969) 
and  homogeneous  distribution  of  blood  and  air  among  the  parabronchi  is 
assumed. 

MacDonald  (1970)  gives  ranges  for  the  diameter  of  the  parabronchial 
lumen  (0.33-0.77  mm)  and  for  the  diameter  of  the  parabronchial  unit 
including  the  lumen  and  capillary  network  (0.77-1.26  ram)  for  the  chicken 
(2  kg,  female).  Payne  (1960)  gives  ranges  for  both  the  length  of  para- 
bronchi (1.0-4.0  cm)  and  the  diameter  of  parabronchial  lumina  (1.39-1.74 
mm),  also  in  chicken  (2.90  kg,  female).   Based  on  these  data,  representa- 
tive  values  have  been  selected  as  i   =  2.42  x  10   cm,  F  =  4.58  x  10 
cm  ,  L  =  1. 7  cm. 


10   mmol-ml   -torr"  ),  and  the  diffusion  coefficients  of  0.  (0.178 

2-1  2     1 

cm  -sec   ,  Handbook  of  Chemistry  and  Physics)  and  CO.  (0.139  cm  -sec   , 

ibid)  in  air,  the  Krogh's  constants  of  diffusion  in  the  gas  phase  are 

determined  to  be  KQ2=6.28xlO"4  and  K   =4.90xl0~4  mmol-min"1-torr"1-cm"1. 
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Kawashiro,  and  Scheid,  1975)  linear  approximations  have  been  taken  as 

802  ■  8.04  x  10"5  and  BQ02   =  27.0  x  10_5  ramol>ml~1-torr~1. 

To  evaluate  the  solutions  given  by  equations  (29)  and  (30),  it  is 

necessary  to  evaluate  p  and  a    (equation  (20)).   Multiplying  the  numerator 

and  denominator  of  both  these  equations  by  the  total  number  of  air 

capillaries  in  a  single  parabronchus  (see  below)  and  the  capillary 

length  allows  p  and  a  to  be  expressed  as 

D  D 

p  -  -*-       ,   and   a  =  -r-|-  . 

BQp*  P 

It  remains  to  find  D  and  A  .   Eased  on  a  carbon  monoxide  test 
P      P 

Piiper,  Pfeifer,  and  Scheid  (1969)  report  the  0.  diffusing  capacity  in 

chicken  (1.4  kg)  to  be  1.0  ml-min  -torr  .   From  this  the  0.  diffusing 

capacity  of  the  model  parabronchus  is  selected  as  D  „  =  6.38  x  10 

puz 

mmol-min  -torr  .  The  CO,  diffusing  capacity  is  taken  to  be  twenty 

-3        -1 
times  as  large:   D   „  =  1.28  x  10   mmol-min   -torr. 

Duncker  (1972)  reports  the  diameter  of  air  capillaries  in  chicken 

to  be  10  urn  and  also  estimates  that  one  half  the  volume  of  the  blood-air 

capillary  network  is  occupied  by  air  capillaries.   Based  on  these  data, 

at  a  radius  equal  to  0.045  cm  the  total  cross  sectional  area  of  all  the 

air  capillaries  in  a  single  parabronchus  is  estimated  to  be  A  =  0.243 

P 
2 
cm  and  the  total  number  of  air  capillaries  in  a  single  parabronchus  is 


All  the  quantities  necessary  to  evaluate  the  solutions  have  been 
assigned  numerical  values. 
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To  understand  the  distribution  of  oxygen  and  carbon  dioxide 
in  the  avian  lung,  a  theoretical  treatment  of  gas  exchange  in  the 
parabronchus  of  the  avian  lung  is  described.  The  model  used  is 
modified  after  Zeuthen  (1942).   Cases  which  treat  diffusion  through 
the  air  capillaries,  diffusion  through  the  parabronchial  lumen,  and 
diffusion  through  both  these  spaces  are  considered.   In  these  cases 
the  0.  dissociation  curve  is  approximated  by  a  straight  line. 
Another  case  has  been  considered  in  which  diffusion  through  the 
air  spaces  is  neglected,  but  the  O2  dissociation  curve  is  more 
accurately  approximated  by  a  polynomial.  The  relationship  of  the 
partial  pressures  of  02  and  C02  within  the  blood  capillaries,  air 
capillaries,  and  parabronchial  lumen  to  parabronchial  blood  flow 
and  ventilation  is  graphically  shown. 

The  results  indicate  that  the  maximum  partial  pressure  variations 
of  0^  and  CC^  along  the  length  of  an  air  capillary  are  less  than 
one  torr  and  five  torr,  respectively.  Also,  the  resistance  to  gas 
exchange  imposed  by  diffusion  across  the  blood-gas  interface  is  many 
times  greater  (about  250  times  greater  for  02;  from  8.5  to  33  times 
greater  for  C02)  than  the  resistance  imposed  by  diffusion  through  the 
air  capillaries.  The  assumption  of  bulk  flow  only  through  the  para- 
bronchial lumen  yields  accurate  results  at  resting  and  high  ventilation 
rates,  but  axial  diffusion  in  the  lumen  must  additionally  be  considered 
at  lower  ventilation  rates.   Before  entering  the  parabronchus,  the 
gas  in  the  pre-parabronchial  region  shows  a  reduction  in  02  partial 
pressure  due  to  forward  diffusion  toward  the  parabronchus  and  an 
increase  in  C02  partial  pressure  due  to  backward  diffusion  from  the 
parabronchus. 


